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Abstract 

The present work deals with the resolution of the Poisson equation in a bounded domain made of a thin and pe¬ 
riodic layer of finite length placed into a homogeneous medium. We provide and justify a high order asymptotic 
expansion which takes into account the boundary layer effect occurring in the vicinity of the periodic layer as well 
as the corner singularities appearing in the neighborhood of the extremities of the layer. Our approach combines 
the method of matched asymptotic expansions and the method of periodic surface homogenization, and a complete 
justification is included in the paper or its appendix. 
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Introduction 

The present work is dedicated to the construction of a high order asymptotic expansion of the solution to a Poisson 
problem posed in a polygonal domain which excludes a set of similar small obstacles equi-spaced along the line 
between two re-entrant corners. The distance between two consecutive obstacles, which appear to be holes in the 
domain, and the diameter of the obstacles are of the same order of magnitude S, which is supposed to be small 
compared to the dimensions of the domain. The presence of this thin periodic layer of holes is responsible for 
the appearance of two different kinds of singular behaviors. Eirst, a highly oscillatory boundary layer appears in 
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the vicinity of the periodic layer. Strongly localized, it decays exponentially fast as the distance to the periodic 
layer increases. Additionaly, since the thin periodic layer has a finite length and ends in corners of the boundary, 
corners singularities come up in the neighborhood of its extremities. The objective of this work is to provide a 
sophisticated asymptotic expansion that takes into account these two types of singular behaviors. 

The boundary layer effect occurring in the vicinity of the periodic layer is well-known. It can be described using a 
two-scale asymptotic expansion (inspired by the periodic homogenization theory) that superposes slowly varying 
macroscopic terms and periodic correctors that have a two-scale behavior: these functions are the combination 
of highly oscillatory and decaying functions (periodic of period 6 with respect to the tangential direction of the 
periodic interface and exponentially decaying with respect to d/S, d denoting the distance to the periodic interface) 
multiplied by slowly varying functions. This boundary layer effect has been widely investigated since the work of 
Sanchez-Palencia [37, 36], Achdou [2, 3] and Artola-Cessenat [5, 6 ]. In particular, high order asymptotics have 
been derived in [4, 27, 13, 9] for the Laplace equation and in [34, 35] for the Helmholtz equation. 

On the other hand, corner singularities appearing when dealing with singularly perturbed boundaries have also 
been widely investigated. Among the numerous examples of such singularly perturbed problems, we can men¬ 
tion the cases of small inclusions (see [29, chapter 2] for the case of one inclusion and [ 8 ] for the case of several 
inclusions), perturbed comers [16], propagation of waves in thin slots [23, 24], the diffraction by wires [14], or 
the mathematical investigation of patched antennas [7]. Again, this effect can be depicted using two-scale asymp¬ 
totic expansion methods that are the method of multiscale expansion (sometimes called compound method) and 
the method of matched asymptotic expansions (see [38, 29, 22]). Following these methods, the solution of the 
perturbed problem may be seen as the superposition of slowly varying macroscopic terms that do not see directly 
the perturbation and microscopic terms that take into account the local perturbation. 

Recently, Vial and co-authors [39, 11] investigated a Poisson problem in a polygonal domain surrounded by a 
thin and homogeneous layer, while Nazarov [31] studied the resolution of a general elliptic problem in a polyg¬ 
onal domain with periodically changing boundary. In their studies they have combined the two different kinds 
of asymptotic expansions mentioned above in order to deal with both corner singularities and the boundary layer 
effect. Based on the multiscale method, the authors of [39, 11] constructed and justified a complete asymptotic 
expansion for the case of the homogeneous layer. For the periodic boundary in [3 1 ] the first terms of the asymptotic 
expansion have been constructed and error estimates have been carried out. This asymptotic expansion relies on 
a sophisticated analysis of solution behavior at infinity for the Poisson problem in an infinite cone with oscillat¬ 
ing boundary with Dirichlet boundary conditions by Nazarov [30], where he published an analysis for Neumann 
boundary conditions in [32]. In the present paper, we are going to extend the work for the homogeneous layer 
and the periodic boundary by constructing explicitely and rigorously justifying asymptotic expansion for the above 
mentioned periodic layer transmission problem to any order (with Neumann boundary conditions on the perfora¬ 
tions of the layer). 


1 Description of the problem and main results 

1.1 Description of the problem 

In this section we are going to define the domain of interest G K^, its limit when (5 —> 0 and the problem 
considered. With the coordinates x = {xi,X 2 ) of let Ob and Ot be the two adjacent rectangular domains 
defined by 


Ob = (-L, L) X {-Hb, 0) , Ot = {-L\ L') x (0, H^) , 

where L' > L, Hb and are positive numbers. We denote by F the common interface of Ob and Ot, i- e., 

f = dn^r]dn^ and r = (-L,L) X {0}. 
and we consider the (non-convex) polygonal domain (see Fig. 1 a) 

O = Ob U Ot U F , 

which has two reentrant comers at Xq = (±L, 0 ) with both an angle of 

Besides, let Ohoie G be a smooth canonical bounded open set (not necessarily connected) strictly included in 
the domain ( 0 ,1) X (-1, 1). Then, let N* := N \ {0} denote the set of positive integers and let (5 be a positive real 
number (that is supposed to be small) such that 


3 


( 1 . 1 ) 







0 -- 


-ffB- 



X-+ 

^0 

0 - 

^0 


Qb 


-Hb- 




-+- 


-+- 


-+- 


-+- 



-L’ -L Q L L' 

(a) The domain Q — fix H fin H T. 


-L' -L 0 L L' 

(b) The domain of interest 


Figure 1: Illustration of the polygonal domain Q and the domain of interest 17“^. 


Now, let n^oie ^ '^hin (periodic) layer consisting of q equi-spaced similar obstacles which can be defined by 
scaling and shifting the canonical obstacle flhoie (see Fig. lb): 

g 

^hole “ U <^{^hole + (^ ~ l)6l} | ■ (1-2) 

e=i 

Here, ei and 62 denote the unit vectors of and 6 is assumed to be smaller than Ht and such that does 

not touch the top or bottom boundaries of 11. Finally, we define our domain of interest as 

= (Hb u Ht u r)\fi^. 

Its boundary dfl^ consists of the boundary of the set of holes F*^ = and F^ \ F^ = dil, the 

boundary of 11. Here and in what follows, we denote by n the outward unit normal vector of dV,^. Note, that in 
the limit (5—^0 the repetition of holes degenerates to the interface F, the domain fl"^ to the domain Ht U Ob and 
its boundary dfl^ to 90 U F. 

The domain O"^ being defined, we can introduce the problem to be considered in this article: Seek solution to 

r -Au^ = /, in O^ , 

< Vu'^-n = 0, onF^ (1.3) 

[ = 0, on F^) , 

where / S L^(O^). It is natural to search for € Hp^ (O"^) where 

Hp^ (O"^) = {u € iT^(0‘^) such that u = 0 on Fd} . (1-4) 

The well-posedness of problem (1.3) in Hp^ (O^) directly follows from Lax-Milgram theorem: 

Proposition 1.1 (Existence, uniqueness and stability). Let f G L^(0‘^). Then, for any ^ > 0 there exists a unique 
solution of problem (1.3) in Hp^ (O^), and with a constant C (independent of 5) it holds 

Mhi(o^) ^ ^ II/IIl=(o^) • (1-5) 

The objective of this paper is to describe the behavior of as <5 tends to 0. For the sake of simplicity, we shall 
assume that / has a compact support in a subset of Hx with distance (jg > 0 to F. Our work relies on a construction 
of an asymptotic expansion of as 6 tends to 0 . 

Remark 1.2. The construction is for simplicity for the specific geometrical setting, where F is a straight line 
ending in two comers of the polygonal boundary 911, where the angles between F and 911 are both ends at angles 
^ or TT, respectively. Nevertheless, the study may be extended to a polygon H of different angles. 

Remark 1.3. It is worth noting that the choice of the boundary condition imposed on the small obstacles F*^ consti¬ 
tuting the periodic layer, here homogeneous Neumann boundary conditions, has a strong impact on the asymptotic 
expansion. A homogeneous Dirchlet condition would yield to a completely different asymptotic expansion (see for 
instance Appendix A in [17], or [12]). 

Remark 1.4. The smoothness o/Hhoie not required for the existence ofu^ G iT^(ll^). The well-posedness result 
remains valid if Llhoie is a Lipschitz domain. However, we use this assumption in the forthcoming analysis (In 
particular in Proposition 4.5 and Section 6). 
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1.2 Ansatz of the asymptotic expansion 

As mentioned in the introduction, due to the periodic layer, it seems not possible to write a simple asymptotic 
expansion valid in the whole domain. We have to take into account both the boundary layer effect in the vicinity 
of r and the additional corner singularities appearing in the neighborhood of the two reentrant corners. To do so, 
we shall distinguish 2 l far field area located ’ far’ from the reentrant corners Xq and two near field zones located in 
the vicinity of them (see Fig. 2). 



Figure 2: Schematic representation of the overlapping subdomains for the asymptotic expansion. The far field area 
{hatched) away from the comers Xq is overlapping the near field area (light grey) in the matching zone (dark grey). 


1.2.1 Far field expansion 

Far from the two corners Xq (hatched area in Fig. 2), we shall see that is the superposition of a macroscopic 
part (that is not oscillatory) and a boundary layer localized in the neighborhood of the thin periodic layer. More 
precisely, we choose the following ansatz: 


«'(x) 


(n,q)GN2 


where x = {xi,X 2 ), and for (n, q) G 


if|a:i|>L, 

X(f)ui^g(x) + n^^g(xi,f) if|a:i|<L. 


Here x : K > (0,1) denotes a smooth cut-off function satisfying 


( 1 . 6 ) 


(1.7) 


X(i) 


1 if|f|>2, 
0 if|f|<l. 


( 1 . 8 ) 


The macroscopic terms ^ are defined in the limit domain Ht U Hb. A priori, they are not continuous across F. 
As for the boundary layer correctors ^{xi, Xi, X 2 ) (also sometimes denoted periodic correctors), and as 
usual in the periodic homogenization theory, there are 1-periodic with respect to the scaled tangential variable Xi. 
Consequently, they are defined in {—L, L) x B, where B is the infinite periodicity cell (see Fig. 3a): 


^ = {(0,1) xM}\Hhoie. (1.9) 

Moreover, the periodic correctors are super-algebraically decaying as the scaled variable X 2 tends to ±00 (they 
decay faster than any power of X 2 ), more precisely, for any {k, i) G N^, 

lim XX.<, = 0. (1.10) 

IA 2 1—>-+oo 

The macroscopic terms as well as the boundary layer corrector terms might have a polynomial dependence with 
respect to In S: there is N (n, q) gN such that 

N{n^q) N{n,q) 

^1? = E (ln<5)*'«n.9.si and = E 

s—0 s—0 

where Un,g,s and n„ ,j ^ do not depend on <5. 
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Remark 1.5. Here and in what follows, although it might be surprising at first glance, we call far field expansion 
the expansion (1.6), i. e., the superposition of the macroscopic terms and the boundary layer correctors. Besides, 
it should also be noted that, for any k G N, we consider 5 s and as different scales as they 

would be different powers of 5. In fact, we shall see that n and q play a different role in the asymptotic procedure. 
Finally, following Remark 1.2, the consideration of the more general case of two angles of measure a, would yield 
to an expansion of the form (1.6) substituting for 5^'^^ (see [11 ]). 


1.2.2 Near field expansions 


In the vicinity of the two comers Xq (light grey areas in Fig. 2), the solution varies rapidly in all directions. 
Therefore, we shall see that 


’w= E « 


¥+9 






for some near held terms g ± 


(n,g)GN2 

dehned in the hxed unbounded domains 


( 1 . 11 ) 


( 1 . 12 ) 


h- = 1C- \[j {flhole +fel} , n+ = IC+\ [j {flhole - fel} 

^GN ^gn* 

shown in Figure 3b and 3c, where K.^ are the unbounded angular domains 

IC^ = {X = (cos6 »±,sine G /±} G 

of angular sectors /+ = (0, and ' = (-f,7r). If the domain flhoie is symmetric with respect to the axis 
Xi — 1/2, then the domain is nothing but the domain 12+ mirrored with respect to the axis Xi = 0. However, 
this is not the case in general. Similarly to the far held terms the near held terms might also have a polynomial 
dependence with respect to In <5, i. e., for all (n, q) G N^, there is N(n, q) G N such that 

N{n.,q) 


where the functions Un^q^±^s do not depend on S. 


1X2] 

xT 


s =0 


Xi 


xE 

Xi 


9 = TT - 


) = 0 


0 = -- 


Stt 

T 


(a) The periodicity cell S. (b) The domain 12 . (c) The domain 12^. 

Figure 3: The periodicity cell B and the normalized domains 


1.2.3 Matching principle 

To link the two different expansions, we assume that they are both valid in two intermediate areas (dark 
shaded in Fig. 2) of the following form: 

= |x = (xi,X 2 ) G Cl,'/s < d(x, Xq) < 2 , 

where d denotes the usual Euclidian distance. The precise dehnition of the matching areas is not important. The 
reader might just keep in mind that they correspond to a neighborhood of the corners Xq of the reentrant corners 
for the far field terms (macroscopic and boundary layer correctors) and to going to +00 for the near held terms 
(expressed in the scaled variables). 

1.3 Far and near field equations 

The ’ansatz’ being assumed, the next objective is to construct the terms g, g and g ± of the far and near 
held expansions (the asymptotic expansions are justihed later, by proving error estimates). This is by far the longest 
part of the work (Sections 1 to 5). The usual starting point of this construction consists in the formal derivation of 
the near held and far held problems, that is to say problems satished by the near and far held terms. 
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1.3.1 Far field equations: macroscopic and boundary layer correctors equations 

Inserting the far field expansion into the initial problem (1.3) and separating the different powers of 6 (the complete 
procedure, based on the separation of the scales, is explained in Appendix A) gives a collection of equations for 
the macroscopic terms and the boundary layer terms: 


Macroscopic equations The macroscopic terms ^ satisfy 

” ^ in fix U tie, 

I 0 otherwise, 

together with homogeneous Dirichlet boundary conditions on Fd 

'u-i,q = 0 on Td- 

Boundary layer corrector equations The boundary layer correctors satisfy 

r-Axn^,,(xi,X)= inB, 

1 g = 0 on aflhole- 


where, for any p e N, 
(-,5 ^ 

^ dx2^ 

p=0 






+ C<i-2(='i.x) + 2£(„ax.nJ,_,(i,,X). (1.16) 


Here, for any sufficiently smooth function v, [A, u] denotes the commutator between A and v, that is to say 

[A, v\u = A{vu) — vAu = Vu • Vm + uAv. 

Moreover, the smooth truncation functions x± ^6 defined by 


X±(^ 2 ) = 


x(A2 ), if±A2 >0, 
0 otherwise, 


(x±(X2) — 1(±X2>0)X(-’^2)), 


and, for p € N, 


(a:i, 0 =^) = 

0X2 /i-s-0± 0X2 


Note that equations (1.13-1.14), posed in the domains Hx and Hb, do not define entirely the macroscopic terms. 
Indeed, we first have to prescribe transmission conditions across the interface T (for instance the jump of their trace 
and the jump of their normal trace across T). This information will appear to be a consequence of the boundary 
layer equations (Section 2). Then, we also have to prescribe the behavior of the macroscopic terms in the vicinity 
of the two corner points Xq. This information will be provided by the matching conditions (Section 5). 

1.3.2 Near field equations 

The near field equations are obtained in a much more direct way. Inserting the near field ansatz (1.11) into the 
Laplace equation (1.3) and separating formally the different powers of <5, it is easily seen that the near field term 
Ut„ satisfy 


-AxU^ g = 

0 

in 


0 

on O/C*, 

= 

0 



1.4 Outlook of the paper 

The remainder of the paper is organized as follows. In Section 2, we investigate the boundary layer problems. We 
derive transmission condition for the macroscopic term ^ up to any order (Proposition 2.4). We also obtain an 
explicit formula for the periodic correctors H^ ^ (see (2.30)). In particular, we shall see that the periodic corrector 
g is completely determined providing that the macroscopic terms ^ are defined for p < q. 

Then, Section 3 is dedicated to the analysis of the far field problems (consisting of the far field equations (1.13) to¬ 
gether with the transmission conditions (2.29a),(2.29b)). We first introduce two families of so-called macroscopic 
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singularities ^ (Proposition 3.5 and Proposition 3.8). These functions are particular solutions of the homoge¬ 
neous Poisson equations (with prescribed jump conditions across the interface T) that blow up in the vicinity of 
the reentrant comers. These two families are then used to derive a quasi-explicit formula for the far field terms 
(Proposition (3.27)). This quasi-explicit formula defines the macroscopic terms ^ up to the prescription of 2n 
constants called q), 1 < m < n. 

Section 4 deals with the resolution of the near field problems (1.18). As done for the macroscopic terms, we de¬ 
fine two families of near field singularities S^, that are particular solutions of the homogeneous Poisson problem 
posed in that blow up at infinity (Proposition 4.8). Based on these near field singularities, we then derive a 
quasi-explicit formula for the near field terms (4.21). Here again, this quasi-explicit formula defines the near field 
terms (7^ ^ up to the prescription of n constants called q±),l<w<n. 

Section 5 is dedicated to the derivation of the matching conditions and the definition of the terms of the asymptotic 
expansions. Based on an asymptotic representation of the far field terms close to the reentrant corners and of the 
near field terms at infinity, we obtain a collection of matching conditions (5.11),(5.13),(5.14) and (5.15) that permit 
to determine the constants ^ for the near fields and the constants q) macroscopic fields. 

As a consequence, all the terms of the asymptotic expansion are then constructed (through an iterative procedure). 

Finally, Section 6 deals with the justification of the asymptotic. 

2 Analysis of the boundary layer problems: transmission conditions 

This section is dedicated to the analysis of the boundary layers problems (1.18). It permits us to derive (necessary) 
transmission conditions for the macroscopic fields ^ across F (Proposition 2.4). For a given n G N, we shall 
propose a recursive procedure to write the jump of the trace and of the normal trace of ^ as linear combinations 
of the mean values of the macroscopic fields of lower order tangential derivatives. This 

procedure is done by induction on the index q and is completely independent of the index n and of the superscript 
5 (of g). That is why we shall omit the index n and the superscript <5 in this section. 

For any sufficiently smooth function u defined in U we denote by 0)]r and (Mjxi, 0))r its jump and 
mean values across F: for |a;i| < L, 

[M(a;i,0 )]p= lim {u{xi,h) — u(xi,—h)), {u{xi,0))r = - lim {u(xi,h) + u{xi,—h)). (2.1) 

h^0+ 2 h^0+ 

Let (/g)qgN be a sequence of functions belonging to L^(flT) H L^(Ob) that are compactly supported in Ht. We 
consider the following sequence of coupled problems (obtained by rewriting (1.13),(1.14),(1.15) and omitting the 
index n): 

-Aug= fq inHTUHB, r -Axn,(xi,X) = G,(xi,X) 

Uq = 0 onFu, \ = 0 

where 

G,(xi,X) = ^ {2{gpiX)){dP^Uq.p{xuO))r + ^ fe(X)] [dP^Uq.,ixuO)]^) 

+ 2d,,dx,Tiq_i + dl^nq_2. (2.3) 

Here, we use 

(5,(X)):=1[A,x++X-](^) , [5p(X)]:=[A,x+-X-](^),peN, (2.4) 

and later also p^(X)) := [A,x^] will be needed. As previously, we impose Hg to be 1-periodic with 

respect to Xi and to be super-algebraically decaying as |A 21 tends to -foo: for any (fc, £) G N^, 

lim n, = 0. (2.5) 

IX2I—>-1-00 ^ 

Note that the right-hand side Gq in (2.3) corresponds to the right-hand side G^ ^ of Problem (1.15) (for a given 
n). The problems for H^, g G N are coupled to the others by the source terms. In difference, the problems for Uq, 
g G N are not complete and their coupling to other problems will be exposed in following. 


in B , 
on (9Hhoie 


( 2 . 2 ) 



The present section is organized as follows; in Section 2.1, we give a standard existence and uniqueness result 
(Proposition 2.2), which shows that under two compatibility conditions the boundary layer problems for Ilq in (2.2) 
have a unique decaying solution. In Section 2.2, we use Proposition 2.2 to derive transmission conditions for the 
first two terms ug and ui (see (2.12)-(2.18)-(2.19)-(2.27)), and we obtain an explicit tensorial representation for 
the associated boundary layer correctors (cf. (2.13)-(2.20)). Finally the approach is extended in Section 2.3 to 
obtain transmission conditions up to any order for the macroscopic fields Ug. 

Remark 2.1. The asymptotic construction described in this section is entirely similar to the construction of a 
multi-scale expansion for an infinite periodic thin layer (without corner singularity). A complete description of 
this case may be found in [37], [3], [1 ], [35] and references therein. 

2.1 Preliminary step : existence result for the boundary layer problem 

In this subsection, we give a standard result of existence for the boundary layer corrector problems for n, in (2.2). 
It will be subsequently used to construct exponentially decaying boundary layer correctors. Let us introduce the 
two weighted Sobolev spaces 


V±(e) = {nGHi„(6),n(0,X2)=n(l,X2), and (nu;±) GHi(S)} 


( 2 . 6 ) 


where the weighting functions wf{Xi,X 2 ) = x(2f2) exp(±i^^). The functions of V~(B) correspond to the 
periodic (w.r.t. Xf) functions of that grow slower than exp(J^^) as X 2 tends to ±00. By contrast, the 

functions of correspond to the periodic functions of Hl^^{B) decaying faster than exp(—as X 2 tends 

to ±c». As a consequence, they are super-algebraically decaying, means they satisfy (2.5) for i < 1. Note also 
thatV+(6) c V-{B). 

Based on this functional framework, we consider the following problem: for given g G (V“ (B))' find 11 G V~ (B) 
such that 

-Axil = g inB, 


dnB = 0 OnOflhole, 

dxM^.X2) = dx,B{l,X2), A 2 GK. 


(2.7) 


Proposition 2.2. 1. Problem (2.7) has a finite dimensional kernel of dimension 2, spanned by the functions 

A/” = Ig and T), where T) is the unique harmonic function ofV~{B) such that there exists Poo G R 

P(Ai, A2) = P(Ai, A2) - X+iX2){X2 + Poo) - X-iX2){X2 - Poo) 
belongs to V'^{B). The constant Poo only depends on the geometry of the periodicity cell B. 

2. If f is orthogonal to P and Af in the Lf (B) sense, meaning that 


p(X)P(X)dX = 0 

3 

g{X)M{X)dX = 0 


(Cv) 

(Cat) 


then, there exists a unique solution 11 G V'^(B). 

3. Conversely, if problems (2.7) admits a solution 11 G V'^(B), then it satisfies the compatibility conditions 

(Cv), (C^). 

For the proof of the previous proposition we refer the reader to [32, Prop. 2.2] and [15, Sec. 5]. General results 
on the elliptic problems in infinite cylinder can be found in [26] (Chapter 5). Note that all these results remain the 
same with a different exponential growth or decay constant in the definition of wf unless it does not exceed an up¬ 
per bound which is determined by the least exponentially decaying or growing functions in the kernel of — A in B. 

Based on the previous proposition, we shall construct X) in C {{—L, L), V^{B)). Transmission conditions 

for the macroscopic terms [itgjr and \dx^Uq\Y will directly follow from the compatibility conditions (C-d), (Ca^) 
applied to Problem (2.2) for Ilg, g G N. It will guarantee that the boundary layer correctors IVq are exponentially 
decaying. Let us give a couple of useful relations, which are easy to obtain by direct calculations (noting that 
g^ = [A, x±]l)> and will be extensively used in the next subsections: 
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Lemma 2.3. The following relations hold: 


/ {goiX))V{X)dX = 0, 

Jb 

[ {go{X))M{X)dX = 0, 
JB 

[ {gi{X))V{X)dX = V^, 
JB 

[ {gi{X))M{X)dX = 0, 


L 

L 

L 

L 


[go{X)]V{X)dX = -2, 
[go{X)]N{X)dX = 0, 
[gi{X)]V{X)dX = 0, 
[gi{X)]N{X)dX = 2. 


( 2 . 8 ) 


2.2 Derivation of the first terms 

We can now turn to the formal computation of the first solutions of the sequence of Problems (2.2). We emphasize 
that the upcoming iterative procedure is formal in the sense that we shall provide necessary transmission conditions 
for the macroscopic terms Uq but we shall not adress the question of their existence in this part (this question will 
be investigated in Section 3). Throughout this section, we assume that the macroscopic terms exist and are smooth 
above and below the interface F. 


2.2.1 Step 0: [uo]r and Flo 

The limit boundary layer term (or periodic corrector) Ho is solution of 


-Axno(xi,X) = Go(xi,X) in .8,^ 

a„Flo = 0 on dflhoie, 


(2.9) 


where Go (xi, X) = 2(po(X))('Uo(a;i, 0))r+5 [(/o(X)] [uo(a;i, 0)]p. Problem (2.9) is a partial differential equation 
with respect to the microscopic variables Xi and X 2 , wherein the macroscopic variable xi plays the role of a 
parameter. For a fixed xx in {—L,L) (considered as a parameter), Go(a;i,-) belongs to {V~{B)y since it is 
compactly supported. Then, in view of Proposition 2.2, there exists an exponentially decaying solution Flo(xi, •) G 
12+ (B) if and only if the two compatibility conditions (Cp, Cj^f) (Prop. 2.2) are satisfied. Thanks to the second line 
of Lemma 2.3, 

[ Go{xi,X)AfiX)dX = 0, (2.10) 

JB 

which means that is always satisfied. Besides, in view of the first line of Lemma 2.3, 

[ Go{xi,X)V{X)dX = -[uo{xx,0)]^. (2.11) 

JB 


As a consequence, we obtain a necessary and sufficient condition for Ho to be exponentially decaying: 


[uo(xi,0)]r =0. (2.12) 

This condition provides a first transmission condition for the limit macroscopic term uq. Under the previous 
condition, Go(xi, X) = 2(go(X))(uo(a;i, 0))r, and, using the linearity of Problem (2.9), we can obtain a tensorial 
representation of Flo € C {(—L, L), 12+(8)), in which macroscopic and microscopic variables are separated: 


no(a:i,X) = (uo(xi,0))rVFo‘(X). (2.13) 

Here the profile function Wo(X) is the unique function of )2+(8) satisfying 


r -AxVFo*(X) 

= 

in 8, 


< dnW^ 

= 0 

on 90hoie, Fq{X) = 2 {go{X)). 

(2.14) 

{ dx,w^io,X2) 

= dx,W^{l,X2), 

X 2 G ffi, 


A direct calculation shows that 

1U*(X) = (1 - 

■X(2^2)). 

(2.15) 


Note that the continuity of Hq with respect to xi is a consequence of the continuity of Gq with respect to xi. 


2.2.2 Step 1: [d^^uf\T, [ui]r, and Hi 

In view of the general sequence of problems (2.2), the second boundary layer (or periodic corrector) Hi satisfies 

r -Axni(xi,X) = Gi(a:i,X) in 8,^ 

1 a„ni = 0 on oUhoie. 
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where, thanks to (2.15) {dx^W^ = 0), 

Gi(a:i,X) = ^ [ 5 i(X)] [aa; 2 Uo(a:i,0)]p + ^ [goO^)] [ui(a^i,0)]p 

+ F^{X){uiixi,0))r + 2(51 (X)) (9,,uo(xi, 0))r. (2.17) 

As for IIo, Problem (2.16) is a partial differential equation with respect to the microscopic variables Xi and 
X 2 , where the macroscopic variable xi plays the role of a parameter. For a fixed xi in (—L, L), Gi(a:i,-) is 
compactly supported in B, and, consequently, belongs to . Then, thanks to Proposition 2.2, there exists 

an exponentially decaying solution ni(xi, •) € V'^{B) if and only if the two compatibility conditions (Cp), (Cat) 
are satisfied. In view of Lemma 2.3, fg (X), [50 (X)], (51 (X)) are orthogonal to M. Then, the second formula of 
the fourth line of Lemma 2.3 gives 


/ Gi(xi,X)AA(X)dX=[9,,uo(xi,0)]p. 

JB 

Therefore, the compatibility condition (C^/) is fulfilled if and only if 

[a 2 : 2 Uo(a:i, 0 )]p = 0 . 

Next, using the first and third lines of Lemma 2.3, we obtain 

[ Giixi,X)V{X)dX = - [ui(cci, 0 )]p + 2P00 (a,,uo(xi, 0 ))r. 
JB 

Therefore, the compatibility condition (Cd) is fulfilled if and only if 

[mi(xi,0)]p = 2X>oo (9a;2Uo(a^i,0))r- 


(2.18) 


(2.19) 


Under the two conditions (2.18)-(2.19), Problem (2.16) has a unique solution Hi in C {{—L,L),V^{B)) (the 
continuity of Hi with respect to xi results from the continuity of Gi with respect to xi). Using (here again) the 
linearity of Problem (2.16), we can write Hi as a tensorial product between profile functions that only depend 
on the microscopic variables Xi and X 2 , and functions that only depend on the macroscopic variable xi (more 
precisely, the latter functions consist of the average traces of the macroscopic terms of order 0 and 1 on T): 

ni(xi,X) = (ui(xi,0))rIU*(X) + (a,2Wo(xi,0))rIUi''(X), (2.20) 

where lUg is defined by (2.14) and lUf G V+ (B) is the unique decaying solution to the following problem: 

-AxVLi"(X) = Fi’'(X) + ^[ 5 o(X)] in B, 

dnWy = 0 onaUhoie, ( 2 - 21 ) 

dx,Wy{0,X2) = dx,Wy{l,X2), X 2 GM, 


( 2 . 22 ) 


where, 

Fi"(X) =2(5i(X))r and Vl= ( Fy{X)V{X)dX = 2V^. 

JB 

It is easily seen that the right-hand side of (2.21) is orthogonal to both Af and T>. A direct computation shows that 

wy{x) = V{X), (2.23) 

the function T> being defined in the first point of Proposition 2.2. 


2.2.3 Step 2: [dx^ui]T ([M 2 ]r and 112 ) 

We can continue the iterative procedure started in the two previous steps as follows. The periodic corrector 112 
satisfies the following equation 


r -Axn2(xi,X) 

i dnB-2 


G2 (xi,X) in,B,^ 

0 on OUhoie- 


(2.24) 


Here, 


G2(a:i,X) = ^ [5o(X)] Hr + ^ [5i(X)] [d^Mv + 

+ Fy{X) + F*(X)52^(uo)r + F^{X) (2.25) 
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(2.26) 


i^Q(X) and i^"(X) are given by (2.14)-(2.22), and, 

F^iX) = -2{g,iX)) + W^{X) , F^iX) = 2 dx, W^iX). 

In formula (2.25), for the sake of concision, we have omitted the dependence on xi of the macroscopic terms. 
To obtain this formula, we have replaced IIq and Hi with their tensorial representations (2.13),(2.20), we have 
substituted —d^^uo{x, 0^) by d^^uo{x, 0^) using the macroscopic equation (2.2 ) (Amq = 0 in the vicinity of T) 
and we have taken into account the jump conditions (2.12),(2.18) for uq. 

For a fixed xi G {—L,L), it is easily verified that G 2 {xi,-) belongs to indeed, the first five terms 

of (2.25) are compactly supported and the last one is exponentially decaying (more precisely, F 2 and F 2 
belong to Then again, the existence of an exponentially decaying corrector 112(0:1, •) G V^{B) results 

from the orthogonality of G2(o;i, •) with M and T). As previously, enforcing the compatibility condition (Cat) 
provides the transmission condition for the jump of the normal trace of mi across T: 

[dx 2 Ui\r =-^2 (uo)r +-1^2 dxi{dx 2 Uo)r, (2.27) 

where 

= - / F:^{X)N{X)dX, A4" = - / F^{X)J\f{X)dX. (2.28) 

JB JB 

Then, enforcing the compatibility condition {Cjg-) provides the jump [u2]r, and the existence of 112 is proved. 
Naturally an explicit expression of [u2]r and a tensorial representation of 112 can be written (see the upcoming 
formulas (2.29a)-(2.30)), but, for the sake of concision, we do not write it here. 

2.3 Transmission conditions up to any order 

We are now in a position to extend the previous approach up to any order. For each g G N, similarly to the first 
steps, our global iterative approach relies on the following procedure: 

1. We compute the right-hand side Gq{xi,X) of the periodic corrector problem (2.2) of order q\ we write 
Gq as a tensorial product between functions that only depend on the microscopic variables Xi and X 2 and 
functions that only depend on the macroscopic variable a;i. More specifically, the latter functions consist of 
the trace and normal trace of the macroscopic terms of order lower than q and their tangential derivatives 
(see (2.9),(2.17),(2.25)). 

2. We compute the normal jump \dx 2 Uq- 1 (xi, 0)]r by enforcing Gq to be orthogonal to Af, i. e. , to satisfy the 
compatibility condition iGx) (s^e (2.12),(2.19)). 

3. We compute the jump [uq{xi, 0)]r by imposing Gq to satisfy the compatibility condition (Cp) ensuring that 
Gq is orthogonal to T> (see (2.18),(2.27)). 

4. We write a tensorial representation of the periodic corrector Ilg introducing at most two new profile functions 
(see formulas (2.13),(2.20)). 

Applying this general scheme, we can prove the following proposition, whose complete proof is postponed in 
Appendix B.l. 


Proposition 2.4. Assume that the macroscopic terms Uq satisfying (2.2) exist. Then, there exists four sequences of 


real constants Afp, Afp, T>p, Dp such that 



Q 

K(a;i,0)]r = ^7?‘aP^(uq_p(a;i,0))r 

+ ^{dx2Uq-p(xi,Q))T , 

(2.29a) 

p^l 

Q 

p^l 

Q 


[dx2Uq{xi,D)]T = ^7Vp‘+i9P+^(u,_p(a;i,0))r 

+ '^My+^dl^{dx2Uq-p(xi,G))T ■ 

(2.29b) 

p^l 

p^l 



In the previous definition.we have used the superscript t (in I?*, N^) to refer to some constants associated with 
tangential derivatives of the average trace of the macroscopic terms. Similarly, the superscript n (in 7?^, N^) is 
used for the constants associated with tangential derivatives of the average of the normal trace of the macroscopic 
terms. 

Remark 2.5. In the proof of Proposition 2.4 (Appendix B.l), we also prove simultaneously that there exist two 
families of decaying profile functions and belonging to V^{B) such that the periodic corrector 11^ G 
C {{—L, L), V^{B)) admits the following representation: 

9 9 

Uq{xi,X) = J2d^xA^,-p{xuO))rW;iX) + J2dPx;"{dx2Uq-p{xi,0))rW;iX). (2.30) 

p—0 p—1 

The definitions of the functions and and of the constants 7?^, Afq, 7?^, Afq are given explicitely in 
(B.4),(B.6),(B.7), and (B.9). 
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We point out that the periodic correctors Ilg do not appear (explicitly) in (2.29b): they have been eliminated. In 
other words, the resolution of macroscopic and boundary layer problems are decoupled and the construction of 11 ^ 
can be made a posteriori. 

3 Analysis of the macroscopic problems (macroscopic singularities) 

Thanks to the previous section (see in particular Proposition 2.4, reminding that the index n and the superscript S 
have been deliberately omitted in the previous section), we can see that if the macroscopic terms ^ (solution to 
(1.13)) exist, they satisfy the following transmission problems: for any (n, q) G N^, 

~ fn,q 
0)] p Qn.qf 

[9x2ui,qixi,0)]^ = 

Ui,q = 0 

where 

9i,q{xi) = XI 

P^l 

P=1 

As previously mentioned, the constants T)^, which only depend on the geometry of the periodicity 

cell B, are defined in (B.6)-(B.9). 

The present section is dedicated to the analysis of Problems (3.1). In Subsection 3.1, we give general results 
of well-posedness for transmission problems: we first introduce a variational framework, then we present an 
alternative functional framework based on weighted Sobolev spaces. In Subsection 3.2, we explain the reason why 
the variational framework is not adapted for the resolution of Problem (3.1) for g = 1 (and higher). This leads 
us to consider singular (extra-variational) macroscopic terms that may blow up in the vicinity of the two corners. 
In Subsection 3.3, we construct several sequences of singular functions that are used in Subsection 3.4 to write a 
general formula for the macroscopic terms (Proposition 3.11). 

3.1 General results of existence for transmission problem 

The problems under consideration can be investigated using the general framework for transmission problems 
posed in polygonal domains developed in [33]. In the present paper, we first recall a classical well-posedness 
result based on a variational form of the problem. Then, based on weighted Sobolev spaces, we describe the 
behavior of the solutions close to the two reentrant corners. 

3.1.1 Variational framework 

Let us introduce the classical Hilbert spaces associated with our problems 

(Ht U Hb) = {u G H^(r 2 x U Hb), s.t. u = 0 on Tpi} , 

which incorporates discontinuous functions over T (see Figure la). Its restrictions to Ht and Hb are denoted by 
Hp^ (^It) and Hp^ (Hb). We denote by Hqq^(F) the restriction of the trace of the function Hp^ (^^t) to F (for a 
complete description of the trace of functions, we refer the reader to [19].), i. e., 

e H'/^(r),s.t.3u e Hp^(Ht) : -c = ^onrj. 

Naturally, the space HgQ^(r) is also the restriction of the trace of the functions of Hp^(Hb) to F. Based on a 
variational formulation, and thanks to the Lax-Milgram lemma, we can prove the following well-posedness result: 

Proposition 3.1. Let f G g G HQQ^(r), and f) G Then, the following problem has a unique solution 

u belonging to Hp^ (Ht U Hb).‘ 


in Ht U Hb , 


onFu, 


fn,q — 


/ ifn = q = 0, 
0 otherwise. 


(3.1a) 


,{xi,0))r , 


p^l 

Q 


’_g_p(xi,0))r + X-^P+i 

p=i 


(3.1b) 

(3.1c) 


Am = f 


in Hx U Hb , 




3.1.2 Weighted Sobolev spaces and asymptotic behaviour 

In the next subsections, we shall study the behavior of the macroscopic terms in the neighborhood of the two 
reentrant corners. It is well-known that the Hilbert spaces (resp. i/™(OT)) are not well-adapted to this 

investigation. By contrast, the weighted Sobolev spaces provide a more convenient functional framework. We 

refer the reader to the Kondrat’ev theory (see [25], [26, Chap. 5 and Chap. 6] for a complete presentation of these 
spaces and their applications). In this part, we introduce the weighted Sobolev spaces associated with our problem 
following the presentation of [26, Chap. 6]. Let us first define the polar coordinates (r^, 9 ^) centered at the vertex 
Xq, i. e., 

xi - {x.q)i = cos{9^), a ;2 - (xq )2 = sin(0=‘=). (3.3) 

Next, we consider the two infinite angular (or conical) domains )C^± centered at Xq of opening ^ 

= {(r=‘= cos6>"‘',r=‘= sin^"*") G > 0,0='= e/=*=} , 7+= (0, ^), /“ = (-^,7r), (3.4) 


and, for i G {0,1, 2}, we define the space V 2 as the closure of \ {0}) with respect to the norm 




a = (ai,a2) G |a| = ai-I-a2- (3.5) 


Then, let 

Xl (x) = (1 - X(2r=^/L)) (3.6) 

be the cut-off function equal to one in the vicinity of Xq and vanishing in the vicinity of x[^ (the support of Xl 
is localized in the neighborhood of the vertice Xq), and let Xl = 1 ~ Xl ~ Xl ■ remind that the truncation 
function x is defined by (1.8). For £ G {0,1, 2}, we introduce the space V 2 





llxi 




Ilx^ 


3('C^+) 


11x2 




< -fcx) 


(3.7) 


equipped with the following norm 


lkllw„(n) — WXl^Wv^ a{K -) + IIXl^IIv;^ (K +) + IIXL^IIff^(a) • (3-8) 

-P ,P Xq -P Xq 

Here, we have used the convention 77° (fl) = Note, that the space 142^^ (H) is independent of the exact 

choice of x and so the truncation functions Xl that 

Vlp, (H) C VIp (H) for any /3' < /? . (3.9) 

In the same way, we also define (resp. V 2 as well as their associated norm || • Wye (rssp. 

II ■ II ^ (resp. Ob) in the definitions (3.7) and (3.8). Finally, for £ G {1, 2}, we introduce 

the space of the trace of the functions in V 2 ^^(Ot) on the interface F. As norm in we take 


ll“ll<//(r) = 


{ll^^lly7„(nT) : ^ y2,p,'v\r = ■ 


When studying the behavior of the far field terms close to the reentrant corners, the set 

r OiTTi 

A = < Am G K, such that Am = m G Z \ {0} 


(3.10) 


(3.11) 


of singular exponents will play a crucial role (see [19, Chap. 1 - 4]). It consists of the real numbers A whose square 
A^ is an eigenvalue of the operator 

. r V{A) = 771 ( 0 , f) n 772(0, f ) c L2(0, f^ L2(0, f) , 

’ 1 u !-)■ Au = —u" . 


Note that the associated eigenvectors are given by 

Wra[t) = sin(Amf), TO G Z \ {0}. 


(3.12) 


The following proposition, which is a standard result in the literature on elliptic problems in angular domains 
(cf. [33, Theorem 3.6 and Corollary 4.4] for the proof), provides an explicit asymptotic representation of the 
solution of the transmission problems in a neighbourhood of the comers (see also [26, Chap. 6] for a complete and 
detailed explanation of the overall approach): 
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Proposition 3.2. Let (3 < Obe a real number such that l — j3 ^ A. Assume that f G V2°/3(^T)nV2°/3(^B) C L'^lLl), 

0 S C and f) G C Then, the unique solution u G Hp^(f2T) H Hp^(r2B) of 

Problem (3.2) admits the following decomposition: 

u= ^ c^ Wq^Q^±{e^) + , (3.13) 

l<q<f(l-/3) 

wherew^ G V2^(flT)f3V2^j^(flB), u’q^o,+ (i) ='ifg(i) cindWq^Q-{f) = Wq{t—^) (where Wq were given in {3>.\2)). 
Moreover, there exists a constant C independent of u such that 

ll''«'''lly2%(aT) + ll^'''llvf ^(sIb) + X! < C” (llflly2%(nT) + Ilflly^^cnB)) • (3-14) 

The expansion (3.13) is nothing but a modal expansion of the solution u is the vicinity of the two comers. Without 
doubt a similar expansion could be obtained using the technique of separation of variables (see [20, Chap. 2]). The 
sum J2i<q<S:(i-p) tUg,o,±(6*^) is an asymptotic expansion for —>■ 0 whose remainder decays 

faster to zero as any term in the sum. Obviously, due the embedding (3.9) asymptotic expansions of higher order 
in are obtained when (3 is decreased (or |/3| increased). 

3.2 The necessary introduction of singular macroscopic terms 

3.2.1 The limit macroscopic term and its behavior in the vicinity of the comers 

The limit macroscopic term Mq q satisfies Problem (3.2) with f = /op = f G L^(Ox) H L^(Ob), 0 = 0 and h = 0. 
In view of Proposition 3.1, there exists a unique solution Uq q belonging to Hq p^ (fix) H Hp^ (Hb)- Indeed, Uq q 
is independent of 6 (it will be denoted by uop) and belongs to Tfp (SI), since its trace does not jump across T. 

The existence and uniqueness of uop being granted, we can investigate its behavior in the neighborhood of the two 
reentrant corners. Since we have assumed that / is compactly supported in fix, / G V^j^(flT) H V'2°^(nB) for any 
/3 G U. Then, in view of Proposition 3.2, uop has the following asymptotic expansion in the vicinity of the two 
corners vertices x^; for any k G N, there exists uo,o,fe G ^2/3 (^t) H V' 2^^(nB) for any /3 > 1 — such that 

k ^ 

“0,0= ^ "" Wm,o,±(ff'^) + ■uo,o,k, (3.15) 

m—1 

where are real constants continuously depending on ||/||Here again, the expansion (3.15) could also 
be obtained using the method of separation of variables. 

3.2.2 A singular problem defining Mq j 

To illustrate the fact that the macroscopic terms of higher orders cannot always be variational (i. e. belonging 
to JI^(Qt) n Tf^(flB)), let us consider the problem satisfied by Uq investigating the regularity of ^ and 
Hq I defined in (3.1b) and (3.1c) (we deliberately omit the term uf q for a while). In view of the asymptotic 
expansion (3.15) of Mop, 

5o.i ^ co,i,±(r=^)"^^^ and Iiq,! ~ 

as tends to zeros. The constants co,i,± and do,i,± can be explicitly determined (but, there is not need to write 
their complete expression). As a consequence, Pq ^ does not belong to iTp(f (T) and Hq i is not in L‘^(T). It follows 
that we are not able to construct Mq i € (Hx U Hb). However, we shall see that it is possible to build a function 
Uq I that blows up as as tends to 0. Since this function is not in U Hb), we say that this 

function is singular. To distinguish from singular functions, we denote functions in if ^(fix U Hb) as regular (so 
not meaning C°°-regular functions). 

Remark 3.3. The previous analysis explains why, contrary to the case of an infinite thin periodic layer (see [34], 
[15]), it is not possible to construct an asymptotic expansion of the form 

u\xi,X2) = X! (u«(x) + n„(xi, X)) , 

where Un G Hq p (fix) H Hp (Ob) and n„ are periodic functions with respect to Xi exponentially decaying as 
X 2 tends to ±00. 

Remark 3.4. Since it is not possible to construct regular macroscopic terms, we shall construct singular ones. 
Nevertheless, the exact solution is not singular. As a consequence the far field expansion (1.6), which contains 
singular terms, can not be valid in the immediate surrounding of the two comers. Here, a near field expansion 
(1.11) has to be introduced, which replace the singular solution behavior towards the comers in their immediate 
neighborhood. 
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3.3 Two families of macroscopic singularities ^ 

In this section, we introduce two families of functions, that are ^ and sZ^ q for the right and left corner, 
that will facilitate the dehnition of the macroscopic terms. The functions are dehned recursively in q for each 
TO G N \ {0}. The following subsection is dedicated to the definition of q, where the functions q, 
g G N \ {0} are defined by induction afterwards. 


3.3.1 Harmonic singularities s^^.o (w G N \ {0}) 

For any positive integer n, the terms q harmonic in ft. It does not imply that they vanish because we allow 
for singular behaviors in the vicinity of the two corners. The present subsection is dedicated to the dehnition of a 
set of harmonic functions that admit singularities in the vicinity of the two corners . The forthcoming analysis 
is done for the right corner xj but a strictly similar approach may be carried out for the left comer. To start with, 
we exhibit a sequence of harmonic functions g that behave like jjj vicinity of xj, and which are 

regular in the vicinity of Xg . 

Proposition 3.5. Let to G N \ {0}. There exists a unique harmonic function g vanishing onFi) of the form 

S-m,0 = {r^)~^W_ra,0,+ {S^)xt + (3.16) 

where g belongs to Hq{Q). 

Proof Remarking that g belongs to Lf(fl), the proof of Proposition 3.5 directly follows from the Lax- 

Milgram lemma. □ 


It is worth noting that g does not depend on the cut off function xt- Besides, it is easily verihed that si 

r2 ’ 


m,0 


belongs to V 2 n for any ^ > 1 + For instance, it belongs to + i (^t) H (11b)- 

Naturally, for to G N \ {0}, we can also prove the existence of a set of functions sZm 0 form 


S_„,g = (r ) ^^W-m,O,-{0 )XL+S-m, 0 ^ S-mP ^ -^0 (H) ■ (3.17) 

As for uo,o, we shall write an explicit asymptotic expansion of sl„ g in the vicinity of the two corners. Applying 
Proposition (3.2) to the function sl^ g (noting that Asl^ g vanishes for < L/2), we can prove the following 

Proposition 3.6. Let to G N \ {0} and fc G N. Then, there exist a function belonging to V"2^^(11t) H 

V 2 for any /3 > 1 — and k real coefficients iff (sl^ 0)' ^ A 9 A such that 


k 


= ir^r^W.^.O,+ {0+) + }2etistm,o)ir^)^^Wq,o,+ ie+)+rtrn,k,+ - 


9=1 


(3.18) 


Analogously, there exist a function ^ _ belonging to V 2 H V 2 ^{LIb) for any /3 > 1 — and k real 

coefficients f“(sl„ 0 ), ^ Q k, such that 


Moreover, for any (3 > 1 


s 


+ 

—m,0 


k 


9=1 


2(fc+l) 

3 


there exists a constant C such that 


(3.19) 


k 

E 

9=1 


(K^(sl^_g)| + \iq (slm,o)l) + \\'^-m,k,+ \\vli,{,Q.-v) + \\'rlm,k,+ \\vlp(^B) 


+ \Vtm,k,-\\vlfQ.^) + \\rtm,k,-\\vlfQ.^) < 




ASIt) + IIS-mpllvA 


T (^b) 


(3.20) 


The formulas (3.18),(3.19) provide asymptotic expansions of q in the neighborhood of x^. Again, despite their 
apparent complexity, they are essentially modal expansions of sl^ g that can be also obtained using the separation 
of variables. Note that the remainder k + orthogonal to the functions Wqfi^+, for q < k: 


//+ 


r + 

— m,fc,+ 


(r+,6»+)w5,o,+ (6»+)d6»+= 0 Wq < k , 


if is small enough (/. e., where Xl — ^)- case, the coefficients q) can be computed as 


e±t.+ 


(sZ o) = 





1’9,0,± 




(3.21) 


Remark 3.7. It is known [26, Chap. 6] that any function v G V'2^^(nT) H V 2 jj{LlB) for /3 > 1 + ^ satisfying 
Av = 0 in Cl is a linear combination of the functions g, 1 < /c < to. 
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3.3.2 The families m S N \ {0}, q gN \ {0} 

In order to construct the macroscopic terms, it is useful to introduce the family of functions {m, q) G (N*)^ 

(remember that N* = N\{0}), corresponding to the ’propagation’ of st^ o (recursively) through the transmissions 
conditions (3. lb),(3.1c): 


Proposition 3.8. For any {m,q) G (N*)^ there exists a unique function ^ G V 2 H V 2 for 

/3 > 1 + ^ + g satisfying 


—As^ = 

c+ — 


[■S^m,o(2;i,0)]r 


[9x2S^m.g(a;i,0)]r = 


0 in fix n flsj 

0 onVo , 

'^'Dldl^{stm,q-p{xi,0))T +^'Dldl;^{da,^stm,q-p{xiMT , (3.22) 


q 


p^i 

q 


'^■f^p+idlt^{st^,q-p{xi,Q))r + ^Afp+idP^{d^^stm,q-p{xi,0))r, 

p—1 p—1 


which admits the following decompositions 

• For any k gN, there exists a function qk + belonging to V 2 pi (Ob) n 1/^2^ pi (Ox) for any /3' > 1 — 
and real constants if {sf^ g_p), 0 < p < q, 1 < n < {k + 1) + such that 


S-m,q = {r^) '^W-m,q,+ id'^An-r+) 

+ E E C(st™,,-p)(r+)^’'-^’tu„.p,+ (0+,lnr+) + (3.23) 

P=0 i<„<(fe+i) + |p 

where Wnfi,+ {d^■, Inr’*') = Wnp,+ {d^) are given in Proposition 3.2 and, for p > 1, Wn,p,+ {d'^, Inr'*') are 
polynomials in Inr^ whose coefficients (functions of 9'^) belong to C°°([0,7r]) H ^]) (here [a, 5] 

denotes the closure of the intervall (a, b)). 

• For any k G N, there exists a function ^ _ belonging to V 2 pi {^t)FV 2 pi {^b) for any (]' > 1— 

and real constants if q-p)’ 0 ^ P ^ 9. 1 < < (fc + 1) + |p such that 

1 

«-m,9 = E E ^f(s-m,q-p)ir~)^''~^Wn,p,-{e-,\nr-) + (3.24) 

p=0 i<„<(fe + i) + |p 

where Irir”) = Wn,o.-(^~) Prop. 3.2) and Wn^p^-{9~,hir~) are polynomials in lnr“ 

whose coefficients (functions of9~) belong to C°°([0, |]) n ff])- 

The proof of Proposition 3.8 is in Appendix C. It is strongly based on the explicit resolution of the Laplace 
equation in so-called infinite conical domains for particular right-hand sides of the form r^(lnr)", A G K, n S N 
(see Section 6.4.2 in [26] for similar results). The proof consists of constructing an explicit lift of the singular part 
of the jump values (3.23) in order to reduce the problem to a variational one (as already done for q). 

Remark 3.9. In the same way, for each m G N* we can define by induction a sequence of functions (sZm q)qeN* 
as follows: sZm. q fbe unique function belonging to V^p (fix) L {^B)for any /3 > 1 -f ^ Fq that satisfies the 
transmission problem obtained from (3.22) by substituting q-p for sZm q-p tbe jump conditions, and the 
asymptotic expansions obtained interchanging (3.23) and (3.24), replacing the superscripts plus by superscripts 
minus. 


3.3.3 Annotations to the singular functions 

Let us comment the results of the previous proposition and of Proposition 3.6: 

- For m > 0 fixed, the family q)geN provides particular singular solutions to (3.1). 

- The exponents A of r~ and r+ appearing in the asymptotic expansions (3.23),(3.24) are singular exponents 

A G A as well as ’shifted’ singular exponents of the form A = A„ —p, A„ G A, where the integer p is between 
1 and q. The most singular part of q in the vicinity of xj is ^ _|_(0+, lnr+), while 

the most singular part of ^ in the vicinity of Xq is (d~,lnr~). Consequently 

q singular’ in the vicinity of the left corner than in the vicinity of the right one. 
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- The function ^ depends only on the functions p for p < q. In others words, providing that n ^ m, 
the definition of the families q,q & N}, {sl„ g G N} are independently defined. 

- The functions Wn,q,± are defined in (C.4),(C.16). However, for the forthcoming derivation of the asymptotic 
expansion and its analysis their explicit expression is not important. Even so these functions will appear 
again in the definition of the near field singularities (see Lemma 4.3). 

- Problem (3.22) alone does not uniquely determine the function Indeed, in view of Remark 3.7 the 

solution of (3.22) is defined only up to a linear combination of Q,n < m + l^}. However, imposing 
additionally the singular behavior close to the corners given by (3.23) and (3.24) (using the fact that the 
functions Wn,q,± are uniquely determined) restores the uniqueness (cf. Remark 3.7). 

- For a given fc G N, the constants 1 < 9 < fc, and the remainders q k ± satisfy an estimate 

of the form (3.20) that has been omitted for the sake of concision. 

- The constants g) are intrinsic to the singularity functions and can be obtained by similar kind of 

formulas as (3.21). Each constant g) appears in the decomposition of several singularity functions 

in (3.23) and (3.24). 


3.4 An explicit expression for the macroscopic terms 

These part is dedicated to the derivation of a quasi-explicit formula of the macroscopic terms ^ by introducing 
particular solutions to Problem (3.1). As mentioned before, we shall allow ^ to be singular. In view of the 
previous construction, we shall impose that 

2n 

'^i,q € n Vlpin-B.) for any/3 > 1 + Y + g, n > 0 , 

ui,q G n for any/3 > ^ + <7 • 


3.4.1 The macroscopic terms u^ g, g G N 

We remind that the limit macroscopic field ugp (remember that Uq g = uo,o G iTg (H)) satisfies Problem (3.2) 
with f = / G n L^((2 b), 0 = 0 and 1) = 0 (see Section 3.2.1). In this subsection we define the functions 

Ug q = ug^q in the large class of possible singular solutions of (3.1) by imitating the iterative procedure of the 
previous subsection for the definition of the singular functions g (i. e., by ’propagating’ ug^ (recursively) 
through the transmission conditions (3.lb),(3.1c)), by which in turn no additional singular functions are added. 

Proposition 3.10. For any g G N* there exists a unique function Ug g = Ug^q G p (fix) H p (Hb ) for /3 > | + g 

o/ (3.1) which admits the following decompositions: 

• For any A: G N, there exists a function ^ _l_ belonging to V 2 pi (Hb) H V 2 pi (fix) for /3' > 1 — and 

real constants {ug^q-p), 1 < p < q, 1 < n < {k + 1) + ^p such that 

Q 

^o,q = Yl H ^n(wo.5-p)(r+)^"”^w„,p,+ (6»+,lnr+) + r+^,,_^(r+,6»+) . (3.25) 

P=0 l<„<(fc + l)+3p 

• Analogously, for any fc G N, there exists a function Vgg ^ _ belonging to (fix) H V 2 pi {VIb) for any 

jd' > 1 — and real constants £f(ug^q-p), 1 < p < g, 1 < n < (fc + 1) -f |p such that 

q 

^o,q = J2 C(wo,9-p)(r")^""^w„,p-(6»-,lnr") + r+^ ,^__(r-,6»-) . (3.26) 

P=0 l<n<(fc+l)+§p 

Note that the functions Wn,p,± in (3.25) and (3.26) were used already in Proposition 3.8 and are defined in (C.4),(C.16). 
Similar to the singular functions the constants £^{ug^p) are intrinsic and fixed, when ug.o is fixed. From now on 
we consider the macroscopic terms Ug q = ug^q to be defined by Proposition 3.10. 

3.4.2 The macroscopic terms n G N \ {0}, g G N 

We construct g as follows: 
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Proposition 3.11. Let n > 0. For any p G N, let p), 1 < k < n be 2n given real constants. Then, the 

family of functions 

q n 

<q = 9 e N, (3.27) 

± p—0k—1 

satisfies the family of macroscopic problems (3.1). Moreover, the term belongs to H V"2^^(r2B) for 

any /? > 1 + ^ + <7. 

We remind that the functions ^ are defined in Proposition 3.5 for g = 0 and Proposition 3.8 for g > 0 and 
sZm q defined in (3.17) {q = 0) and Remark 3.9 {q > 0). 

Proof The ’most singular’ part of ^ n > 0, defined by (3.27) corresponds to o)'®-n which, in 

view of Proposition 3.8 belongs to V2^^(riT) H for any /3>l + ^+( 7. Asa consequence ^ belongs 

to n V 2 fj{LlB) for any /? > 1 + ^ + q. Next, let us show by induction on q that the family q)geN 

is a particular solution to the family of problems (3.1). The base step (q = 0) is trivial. For the induction step, it 
is clear that ^ is harmonic in fix and in IIb and fulfills homogeneous Dirichlet boundary conditions on Td. It 
remains to show the jump conditions across T. Substituting the definition (3.27) into (3.1b) we can assert that 



± p=0 k—1 


{^—k,q—p—r}r ^xi {^X 2 ^ — k,q—p—r}r') 


r—1 


Interchanging the sum over r and p, using the induction hypothesis, we get the expected jump: 


q—r n 


r—1 ± p—0 k—1 


)(s_fe,,_p_.)r = E d:, + T^r d:p) Kq-r) 


r—1 


r • 


The condition for the normal jump follows accordingly, and the proof is complete. □ 

Let n > 0 and g € N be fixed. The function g (defined by (3.27)) is determined up to the specification of the 2n 
constants q), \ < k < n (There are 2n degrees of freedom). The matching procedure will provide a way 

to choose these constants in order to ensure the matching of far and near field expansions in the matching areas. 


3.4.3 Expression for the boundary layers correctors 

Assume now that g is defined by (3.27).Then, inserting this definition into the formula (2.30) defining the 
boundary layer correctors 11^ g, we find them to be given by 

= Ho., = E %‘'9Si(^o.,-p(a;i,0))r + E wo.g-p(a;i, 0))r , 

p=0 p—1 

and, for n > 0, 


n 


s 

n,q 


n / q q-p 

± k—1 \p—0 r—0 


+ iKT.^-ki<r)dV 

p—1 r—0 


{dx2st 


k,q—p—r 



3.4.4 Asymptotic of the far field terms close to the corners 

Thanks to the previous formulas, we have a complete asymptotic expansion describing the behavior of both macro¬ 
scopic and boundary layer correctors terms in the vicinity of the reentrant corners: for any k G N there exists a 
function Un,q,k,+ belonging to V’2^^(nT) H V2^^(nB) for any /3 > 1 — such that 


+ {0'^, Inr'*") + Un,q,k,+ , 

r=0 m=—n 

where, for any (n, m) € x Z, 

_ n j 

^ ^ ^ ^ fc (^n,p)^m —fcj—p) ? Ti > 0. 

± fc—max(l, —m) p—0 


(3.28) 


(3.29) 


Here, we have used the convention that wo.r.i = 0 for any r € N, that f p) = 0 for any m < —k and p e N. 

Moreover, the notation Yl!m=-n denotes the sum over the integers m G Z such that —n<m<k+p (the integer 
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index m does not exceed [k + |rj, where [aj and [a] denote the largest integer not greater or the smallest integer 
not less than a, respectively). 


For TO < 0, the expression of an,j,m,+ can be simplified since (s_^ q-p-r) vanishes, and q-p-r) — ^ 

unless k = —m and q — p — r = 0 (in the latter case, o) = !)• 

(3.30) 

One can also give an asymptotic expansion for 11^ ^ for (xq)i — xi sufficiently small (we remind that Xq = 
((xq)i, (xq) 2) denotes the coordinates of the vertice Xq): 


q fc+|r ^ 

^n,q ~ y ^ ^ ((^o)l ~ ^ ^n,q-r,m,+ Pm,r,+ ~ ^) 3" ^n,q,k,+ ^ 

r—0 m——n 


(3.31) 


where. 


r r 

y, y) = XI ffm,r-p.p.+ (lnf) (y, y) + X 5m.r-p.p,+ (lni) (y, y) ■ (3.32) 


p=0 


P=1 


The functions g^r q + glnr q + polynomials in Inf. Their definitions are given in (C.10),(C.l 1). The 
remainder Tin.q,k,+ can be written as 


p=0 P—0 


(3.33) 


where one can verify (using a weighted elliptic regularity argument, see [26, Corollary 6.3.3]) that the functions 
^n.g.p and belong to Vlp{^T) C for any /3 > 1 - 


Naturally, similar asymptotic expansions occur in the vicinity of the left comer. 


4 Analysis of the near field equations and near field singularities 

The near field terms g ± satisfy Laplace problems (see (1.18)) posed in the unbounded domain 17^ (defined 
in (1.12)) of the form 

{ —Art = / in 17^, 

u= 0 on 9/C"*", (4.1) 

dnu = g on f = 9f7"^ \ 9/C"^. 

In this section, we first present a functional framework to solve the model problem (4.1) (Subsection 4.1). We 
pay particular attention to the asymptotic behavior of the solutions at infinity (Proposition 4.5). Based on this 
result, we construct two families S^, q G N*, of ’near field’ singularities, /. e., solutions to (4.1) with / = 0 but 
growing at infinity as (/?"*") ^ (Subsection 4.2). Finally, we use these singularities to write a quasi-explicit formula 
(see (4.21)) for the near fields terms g ± (Subsection 4.3). Here again, most of the results are explained for the 
problems posed in f7+ but similar results hold for 17“. 


4.1 General results of existence and asymptotics of the solution at infinity 

4.1.1 Variational framework 

As fully described in Section 3.3 in [11], the standard space to solve Problem (4.1) is 


53(17+) = {ve iTioc(17+), Vu e L^(17+) 


l + (i?+)" 


e L^(17+),r; = 0on9/C+ V , (4.2) 


which, equipped with the norm 


laj(n+) 


VT+Wf 


1/2 






(4.3) 
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is a Hilbert space. The variational problem associated with Problem (4.1) is the following: 

find uG 03(11+) such that a(u,t;) = / f{X)v{X)dX+ [ g{X)v{X)da, Vu G 03(H+), 

where o(u, v) = Vu(X) ■ Vv(X)dX. It is proved in [11, Proposition 3.6] (cf. also [39, Lemma 2.2]), that 
f giX)v{X)da< + 

"'Tie 

and that the bilinear form a is coercive on 03(11+) (the seminorm of the gradient |! Vt;|| is a norm on 03(H+)). 

As a consequence, the following well-posedness result holds: 

Proposition 4.1. Assume that 1 + (i?+)^/ G L^(H+) and (1 + (i?+)^)^/^p G i^(rj|(,|g). Then, Problem (4.1) 
has a unique solution u G 03(H+). 

4.1.2 Asymptotic expansion at infinity 

As usual when dealing with matched asymptotic expansions, it is important (for the matching procedure) to be 
able to write an asymptotic expansion of the near fields as i?+ tends to infinity. In the present case, because of 
the presence of the thin layer of periodic holes this is far from being trivial: there is no separation of variables. 
However, Theorem 4.1 in [32] helps to answer this difficult question. 

For the statement of the next results, we need to consider a new family of weighted Sobolev spaces. For f G N (in 
the sequel, we shall only consider £ G {0,1,2}), we introduce the space 03^ .^(^2+) defined as the completion of 

67)?° (fl+) with respect to the norm 

= E 11(1 + P = 1 + (1 + i?+)|0+ - ttI. (4.4) 

p—0 

The norm |j • |j2j2 (fj+j is a non-uniform weighted norm. The weight varies with the angle 0+. Away from the 

periodic layer, i. e., for |0+ — 7r| > e for some e > 0 and i?+ sufficiently large, we recover the classical weighted 
Sobolev norm y^(/C+) (cf. (3.5)): 


lkllvg^(K:+) — 


1/2 


ElK^^) 


t-£+p 




L2(/C+) 


\p—0 


(4.5) 


Indeed, in this part p ~ 1 -|- i?+ for i?+ —> oo. In contrast, close to the layer, i. e., for 0+ —> tt for i?+ fixed, we 
have p —> 1, and the global weight in (4.4) becomes (1 -|- R+'jP-^-^p.o^ 

In the classical weighted Sobolev norm (4.5), the weight (i?+)^ depends on the derivative (p = 0 or p = 1) 
under consideration. It increases by one at each derivative. This is linked to the fact that the gradient of a function 
of the form (i?+)^p(6*), which is given by (i?+)^ ^ {Xg{9)er + g'{9)eg), decays more rapidly than the function 
itself as i?+ tends to +oo (comparing (i?+)^ ^ and (i?+)^). This property does not hold anymore for a function 
of the form (Ar+)^p(26+, 26^) where X+ = (26+,X^) = i?+(cos0+,sin0+) and g G V+(;B) (g is periodic 
with respect to A+ and exponentially decaying with respect to X^)- Indeed, in this case 

V ((X+)^p) = (A(X+)^-ip + (A+)^9^+p) ei + (A+)^9^+p 62, 

which does not decrease as (26+)^“^. This remark gives a first intuition of the necessity to introduce a weighted 
space with a weight adjusted in the vicinity of the periodic layer, i. e., for 0+ —>^ tt. Note that in the case of Dirichlet 
boundary conditions on the holes, the appropriate weighted space to consider is slightly different (see [30]). 


To be used later we prove the following properties of these new function spaces. 

Proposition 4.2. Let 7 G (|, l), p G {1, 2}, A G K, G N. Let x and x± be the cut-ojffunctions defined in (1.8) 
and (1.17). 

- The function Vi = x^^'^ (X^) X- ("^i”) (1’^ i?+)‘^ belongs to 03^ providing that 

/3<x-X + l/2. 
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- The function V 2 = belongs to providing that 

/3 < 1 - A. 


- Let w = w{X^ ^ X 2 ) be a 1-periodic function with respect to X^ such that < + 00 . 

Then, the function V 3 = |Arj^|)'^w(AC^, AT^) belongs to providing that 

p < 7 - A + 3/2. 


- Let w G V^{B) n Hf^yB) such that 

J + |i9^+wp + < + 00 . 

Then, the function V 4 = X-(X;j^)|Xj^|'''“^(lni?+)‘^u'(Xj^, AT^) belongs to providing that 

;3<7-A + 1/2. 

In absence of the periodic layer the solutions of the near field equations might be written as linear combina¬ 
tion of harmonic functions Wm,o,+ (ln 0+), m G N* for i?+ —00 where Wm,r,+ have been de¬ 

fined (C.4). With the periodic layer the behavior far above the layer remains the same, but has to be corrected by 
-^ 1^1 ) (with Pm,r,+ defined in (3.32)) to fulfill the homogeneous boundary conditions 

on Thole- This correction is not harmonic and has a particular decay rate for i?+ —>^ 00 . It can be (macroscopically) 
corrected by and in the neighbourhood of the layer by \X^y'^~^pm,i,+ i}-V\\X^\, X^ X 2 ). 

Then, through a consecutive correction in the form (the cut-off function x has been defined in (1.8) 

^ {{R+f--^Wra,r,+ {\rvR+,e+)x{Xt)+x{X+)\X+y--'^Prr,,,,+ {\w\XtlX+,X+)) (4.6) 

r—0 

the Laplacian becomes more and more decaying for i?+ —> oo, where any decay rate can be achieved, which 
becomes, at least formally, zero forp — oo. The previous observation will be justified in a more rigorous form in 
the following lemma which turn out to be very useful in the sequel. 

For this let us introduce a smooth cut-off function Xmacro,-i- (see Fig. 4) that satisfies 

ixiXt) forX+<-l, 

Xmacro.+ (X+, ^ 2 +) = <^ 1 for X+ > - A , (4.7) 

[1 forX+>-l,|X2+|>3, 

and for m G Z \ {0} the asymptotic block (we adopt this notion from [32]) 

P . 

Um,p,+ {x+,xy) = x(i?+)^ (xmacro.+ (X+,X 2 +) (7?+)^™-’-(i?+, 0+) + 

X-iXt)\Xty--^Pm,r.+ MXtlXt,X+)'^ , (4.8) 
where the cut-off function x- has been defined in (1.17). 

Lemma 4.3. Under the condition 7 G (1/2,1) the Laplacian AUm,p,+ of the asymptotic block Um,p,+ belongs to 
23^ for any P that satisfies 

P<2-Xm+p. (4.9) 

The functions Pm,q,+ are polynomials in In |X/" |. They are periodic with respect to X^ and exponentially decaying 
as Xy tends to ±c». The functions Wm,q,+ are polynomials in In i?+. Note that, for i?+ large enough, if AT/" < —1, 

p 

Ura,p,+ = {{R'^)^'--^Wm,r,+ pnR+,e+)x{Xt) + \X+Pm,r,+ MX+X+, X+)) (4.10) 

r—0 

and, if > — 1, 

p 

Um,p,+ = Y.{R'^y--'^Wra,r,+ ■ (4.11) 

r—0 

The usage of the cut-off functions Xmacro,-i-. X-{Xi), and x{R^) in (4.8) is simply a technical way to construct 
function defined on the whole domain 12+ (and more precisely belonging to i7[Q;,(r2“*’)) that coincides with Um,,p,+ 
for large i?+ as given in (4.10) and (4.11). 
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Figure 4: Schematic representation of the cut-off function Xmacro,+ defined in (4.7) 

Remark 4.4. In view of Proposition 4.2, reminding that ^ ^t) is proportional to (1 —x(X^)), 

for any j G (1/2,1), the asymptotic block lAm,p,+ belongs tofor any fj < 1 — Am- 

Defining Xmaao,-{Xf, Xf) = Xmacro,+ {—Xf,Xf), we Can also define the asymptotic blocks associated with 
as follows; 

P , 

X+iXf)\Xf\>'---p^,r,-MXf\,Xf,Xf)y (4.12) 

We are now in a position to write the main result of this subsection, which proves that for i?+ large and for 
sufficiently decaying right-hand sides, the solutions of Problem (4.1) can be decomposed into a sum of radial 
contributions corrected by periodic exponentially decaying correctors in the vicinity of the layer of equispaced 
holes. In the following, a real number /3 is said to be admissible if /3 — 1 ^ A. 

Proposition 4.5. Let fc G N and 7 G (1/2,1). Assume that f G .^(r2“*')/or some admissible (3 > max(3,1 + 
2(fc+i) ^ + {R '^)^/ € A^(n'^)) and that g G iF^'^^(rhoie) A compactly supported. Then, the unique 

solution u of Problem (4.1) belongs to QJ(n^) H ..^,(12''’) for any (3' < A. Moreover, u admits, for 7 — 1 
sufficiently small, the decomposition 

2 

u = ^.^-n{u)U_ri,p(n),+ + U, p(n) = max(l, [-(/c - n)]), (4.13) 

n—1 

where the asymptotic blocks hl-n,p{n),+ defined in (4.8), 1 < n < k, denote k constants, and, the 

remainder u G for any (3^ such that 13^ < (3. In addition, there exists a constant (7 > 0 such that 

k 

^ ,^(n+) + ll5llrri/2(rM,)) • (4.14) 

n—1 

We remind that [a] (used in (4.13)) stands for the smallest integer not less than a. 

The proof of Proposition 4.5, postponed in Appendix D.3, deeply relies on successive applications of the following 
lemma, which is a direct adaptation of Theorem 4.1 in [32]. To long to be presented in this paper, its proof requires 
the use of involved tools of complex analysis that are fully described in [32]. 

Lemma 4.6. Let 7 G (|, 1)- Let (3^ and 13“^ be two admissible exponents such that (3^ < /3^ and /3^ — /3^ < 1. 
Assume that u G ..^(12’'’) satisfies Problem (4.1) with f G ^^^2 ond a compactly supported g G 

2T^^^(rhoie)- Then, for 7 — 1 < 0 sufficiently small u admits the decomposition 

u= ^ CkUk,3,+ +u, (4.15) 

§(l-/32)<fc<f(l-/3i) 

where the asymptotic blocks Uk, 3 ,+ are defined in (4.8) and for any admissible j3^ G {j3^,f3'^) the remainder 
u G 23|o .y(A2'''). In this case, there exists a positive constant C > 0 such that 

(n+) ^ — C'M|m||qj 2 ^ (n+) + II/IIqjOj (n+) 

§(l-/ 32 )<fc<f(l-/ 3 i) V ^ fir, / 
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We emphasize that the powers of i?+ (or X^) appearing in (4.13) (see (4.8)) are of the form —^n — q, n G N*, 
g S N. Thus, they coincide with the ones obtained for the far field part (see for instance Proposition 3.8). More¬ 
over, as can be expected, the ’leading’ singular exponents ^ correspond to those of the problem without periodic 
layer (here again, as for the macroscopic terms). 


Remark 4.7. The assumption on g in Proposition 4.5 and Lemma 4.6 could be weakened by using the trace spaces 
associated with the weighted Sobolev spaces (see [32]). 


4.2 Two families of near field singularities 

This subsection is dedicated to the construction of two families of functions m G N*, hereinafter referred to 
as the near field singularities for the right and left corner, satisfying the homogeneous Poisson problems 


f-A5± = 0 inf]±, 

J = 0 ond/C^, (4.17) 

[ d„s^ = 0 on \ dJC^ , 

and behaving like for large, where Wm,o,+ = sin(^0+) and Wm,o,- = sin(^(0“ — f)) 

(defined in Proposition 3.2). 

Proposition 4.8. There exists a unique function G .^(LR)for any /3 < 1 — Am and 7 G (1/2,1), satisfying 

the homogeneous equation (4.17) such that the function 

S+=S+-Um,ii+xR,+ (4.18) 


belongs to Q3(n“*’). Moreover, for any k G N*, choosing 1 — 7 sufficiently small, there exists a function TZm.k G 
5j|o j(^^)for any admissible /3° < 1 + and k constants (1 < fc < n) such that Sif admits the 

decomposition 

k 


~ ,+ 'y ^ 2(fc-n) ^ ^ + 'R-m,k- 

n—1 

In addition, for any /? < 1 — Am, there is a constant C > 0 such that 


(4.19) 


k 

+ \\'^m,k\\^2 ^g(f2+) < C' ||5'm||2j2^^(Q+) • 

n—l 


(4.20) 


In the same way as the near field singularities for the right corner 5+, the near field singularities for the left corner 
S'” can be defined. Note that for m > 2 there are several functions satisfying the homogeneous equations (4.17) 
and behaving like {R'^)^^Wm,o,+ (leading term) for large i?+. Indeed, admitting the existence of the functions 
Sm, any function of the form 5+ + Yl'k=i would also fulfill these requirements. Nevertheless, the (4.19) 

restores the uniqueness by fixing (arbitrary) Ok, k = 1 ,... ,m — 1 to 0. 

Proof. The proof is classical and is very similar to the proof of Propositions 3.5, 3.6 and 3.8. We first prove the 
existence of 5"+. The function satisfies 

r -A5+ = /m in f]+, 

< 5+ = 0 ona^+, /m =-A^7m,^l+A„,],-r• 

[ = 0 ondh+\dJC+, 

In view of Lemma 4.3, /m belongs to 9J2 /3(^^) for /3 < 2 — Am + [1 + Ami- Noting that 2 — Am + [1 + 
Ami Proposition 4.1 ensures the existence and uniqueness of 5"+ G 93(n“''), and, hence, the existence of 5+. 
Uniqueness of follows directly from the fact that difference of two possible solutions is in the variational space 
93(U+) and satisfies (4.17). Finally, the asymptotic behavior for large i?+ results from a direct application of 
Proposition 4.5 to the function — Um,p,+ choosing p sufficiently large so that A(,5m — blm,p,+) belongs to 
7(11^) for a real number /3 > max(3,1 + (which is, thanks to Lemma 4.3, always possible). □ 


4.3 An explicit expression for the near field terms 

As done for the macroscopic terms in Section 3.4, we can write a quasi-explicit formula for the near field terms 
j-. We shall impose that the functions ^ do not blow up faster than for i?+ —> 00. Since ^ 
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satisfies the near field equations (1.18), it is natural to construct U^g ± as a linear combination of the near field 
singularities S^, 1 < k < n, namely 

n 

Ulg^^=Y,^k{Ulg^±)St, ( 4 . 21 ) 

k=l 

where q ±) are constants that will be determined by the matching procedure and that might depend on <5. 

Naturally, the functions Un,q,± (defined by (4.21)) satisfy the near field equations (1.18) and belong to 
for any /? < 1 — A„ and 7 € (1/2,1). It is worth noting that the definition (4.21) implies that 

K<i,± =0 e N. (4.22) 


Asymptotic behavior for large To conclude this section, we slightly anticipate the upcoming matching 
procedure by writing the behavior of g ^ at infinity. Thanks to the asymptotic behavior of for —> 00 

(Proposition 4.8), we see that, for any K G W, 


U, 


n K 

n.,.+ = E E 
k—1 l——k 


^k{Ut g +)-^-l{S^)U ; |- 2(iC-i) ^^ ^ + Tln,q,K,+ , 


(4.23) 


where TZn.q,K,+ € ^2 ,go(^^) /3° < 1 + Here, for the sake of concision, we have posed 

) = 0 for any r € N such that 0 < r < g — 1, 

and = 0 for any j € N. Then, substituting the asymptotic blocks U-i^k+ 2 -i,+ for their explicit expres¬ 

sion (4.8), we obtain the decomposition 

K 


t/E+=Xmacro.+ (X) ^ ^ (i?+) ^ (i?+, 0+) 


l ——n r=0 

r 2(iC-0-| 


K 

+ X-{X+)Y^ E <,,-/.+ |X+|-T->_,,,,+ (ln|A+|,A+,A2+)+7^„,,.K.+ 

l ——n r—0 


(4.24) 


where we have used the convention tuo,r,+ = 0 and po,r,+ = 0 for any r G Z. Here, 

n 

V(n,<7,/)GN2xZ, Alg,i,+ = E ^k{Ulg^+)J^iiS+). (4.25) 

/t=max(l,Z) 

Note that for n = 0, ^ ; _|_ = 0. Moreover, for I > 0, ^i{Sk) = 5i,k, and consequently 

<,,/.+«,.+)• (4.26) 

Finally, with the change of index —l—^m and summing up over n and q, we can formally obtain an asymptotic 
series of the near field; For < — 1, 


(n,g)GN2 

n 

(n,(3)eN^ m= —oo reN 

(4.27) 


and, for > — 1, 


E ' 

(n,g')GN2 




= S 

(n,g')GN2 


-+q 


E 

m— — oo 




^ (i?+) ^(0+, In i?+) 


(4.28) 


reN 


The near field terms U^g _ can be decomposed in strictly similar way by substituting formally the superscript 
plus into a superscript minus in (4.24) and (4.25). Here, it should be noted that the way to compute 5 - ) is 

different how ^ _i_) are computed. 
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5 Matching procedure and construction of the far and near field terms 


We are now in the position to write the matching conditions that account for the asymptotic coincidence of the far 
field expansion with the near field expansion in the matching areas. Based on the matching conditions, we provide 
an iterative algorithm to define all the terms of the far and near field expansion (to any order), which have not been 
fixed yet. 


5.1 Far field expansion expressed in the microscopic variables 

We start with writing the formal expansion of the far field 6^3+<i{ui^^{x)x{x2/6) + ni gixi, f)) 

(cf. (1.7)) in the matching area located in the vicinity of the right corner (i. e. for small r+). Collecting (3.28) and 
(3.31), summing over the pair of indices {n,q) € and applying the change of scale x+/(5 = (x —Xq)/() X+ 
and so r+/(5 —> we formally obtain for Xi < —1, 


(wl5(x)x(y) + |)) = 

(".9)6N „ 

^ ^ aL„,,,„,+ ^((i?+)^-’-tt;™,.,+ (0+,lni?++ln5)x(X2+) 

(n,(3)eN2 m=-oo reN , 

+ \X+\"-^--pm,r,+ MX+\+\n6,X+,X+)j , (5.1) 

and for Xi > — 1 

n 

(n,q')eN (n,g)eN 2 m=-oo reN 

Note, that the coefficients ^ defined in (3.29), depend for n > 0 on (5 only through the constants p), 

which we are going to fix in the matching process. In the equations (5.1) and (5.2) the terms Wm,r,+ and Pm,r,+ 
appear with a second shifted argument, i. e.. In 7?+ + ln(5 instead of lni?+ and In + In 5 instead of In |2f 
The following lemma is a reformulation of these terms as linear combinations of non-shifted ones and will prove 
very useful in the matching procedure. It is based essentially on the fact that the terms Wm,r,+ are polynomials in 
the second argument and Pm,r,+ in the first. The proof of the lemma finds itself in Appendix E. 

Lemma 5.1. The equalities 


,, j-(0^,lni?^ + In (5) 


Lr-/2J k 

Y, Wm-Zk,r-2k,±{S^MR^) Y (In (5) , 

k—0 2=0 


(5.3) 


and 


p„,,,±(ln|X±|+ln<),X±) 


Lr-/2J k 

Y Pm-3fc.r-2fc.±(ln|X±|,X±) Y^rn, 2 k,i,± (In <5)* 


(5.4) 


k^O 2=0 

hold true, where for any integer i^k such that 0 < i < k and using the notation Wm, 2 k,± (6l±,lni?±) = 
(In R^yWm,2k,i,± {0'^), the constants Cm, 2 k,i± ore given by 


A r ‘\'TT TT 

Cm,2k,t,± = ly- j Wra,2k,i,±{0^)Wm-3kfl,±{d^)d9^, /+ = (0, —), I~ = (--,7r). 

OTT Ji± Z Z 


(5.5) 


Remark 5.2. With the convention that wo,o,± = 0 it follows that Cm, 2 k,i,± = Ofor k = m/3. Moreover, in view 
of the orthogonality conditions (C.15),(C.23), Cm, 2 k,o,± always vanishes ifk f 0 and Cm,o,o,± = 1- 

Inserting (5.3) into (5.2) and noting that |(m — 3k) — {r — 2k) = |m — r, we obtain 


(n,g)eN2 (n,g)eN2 m=-oo 

L^/2J _ k 

r^N fc—0 2=0 

Then, the changes of indices r — 2k —t r and m — 3k —t m give 

n 

Y E E Xt>0, (5.6) 

(n,q')eN2 (n,g)eN2 m=-oo reN 
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where 


Ln/3J k 

y{n,q,l)Gf^ X Z, ^n,q,m,+ ~ ^ 3fc,q,m+3fc.+ ^ ' C'm+ 3 fc, 2 fc,i,+(l^ J) 

fc—0 i—0 

In particular, for m < 0, thanks to (3.30) (and using Remark 5.2), we have 

L"/3J k 

^n,q,m,+ ~ -^mi'^n,q) ^ ' ^n-3k.q.m+3k.+ ^ ' C'm+ 3 fc, 2 fc,z,+(ln (5) . 

fc—1 i—0 

Analogously, for < — 1, we obtain, 

^^^‘^i'^i,qi^)xix 2 /S) +Uig{xi,x/6) = 

(n,g)GN 

n 

^ 5i"+« d^„,,,„,+ ^((i?+)^-’'tt;^.,,+ (0+,lni?+)x(X2) + |X+|^->™,,.+ (ln|X+|,X+)). 

(n,q)eN^ m= — oo reN 

(5.9) 

The previous two expressions have to be compared with formula (4.27) and(4.28), in which the coefficients 
are still not determined, since the constants q ^), m = l,...,nare not fixed yet. We aim 

to match the expansions in the matching zone and, hence, define these constants uniquely. 


(5.7) 


(5.8) 


5.2 Derivation of the matching conditions 

Arrived at this point, the derivation of the matching conditions is straight-forward. It suffices to identify formally 
all terms of the expansions (5.9) and (5.6) of the far field with all terms of the expansions (4.27) and (4.28) for the 
near field. We end up with the following set of conditions: 

K,q,m,+ = ai_rn,q,m,+ , V(n,g)eN^ and TO G Z, m < 71, (5.10) 

where Af^ q m + q m + Were defined in (4.25) and (5.7). As the coefficients Af^ q m + linear combi¬ 

nations of the constants AfruiU^ q +), to = 1, ..., n and the coefficients ^ ^ _i_ are linear combinations of the 
constants p), to = 1,..., n for some p G N we aim now to obtain conditions between those constants. 

Here, we will proceed separately for the the cases to > 0, to = 0 and to < 0. 


For 0 < TO < 77, using the equality (4.26), we have for any (n, q) gW x N, 

^m{Ui^q,+ ) = ai.m,q,m,+ ■ (5.11) 

For TO = 0, for any {n, q) G N, there is nothing to be matched. Indeed, both left and right-hand sides of (5.10) 
vanish (dn,q,o,+ = 0 because C 3 k, 2 k,i,+ = 0, see Remark 5.2). 


For TO < 0, in view of (5.8) and substituting ^ ^ + for its definition (4.25), we have for any in, q) G N^, 


L(n-m)/3j 


= - E 


a 


n—7n—3k,q,m 


+3fc,+ E Cm+3k,2k,i,+ {^kl6y + E (5.12) 




z=0 




which may also be red as follows: for any (n, m) € such that, n > 1, and 1 < m < n, 


Ln/3J 


=-Y. 


n—3fc,g, —m+3A; 


,+ y]C_™+3fc.2fc.*.+ (ln<5)* + y] (5.13) 




2 = 0 


fe=l 


Here again, we can write similar matching conditions for the matching area located close to the left corner. These 
conditions link the macroscopic terms ^ to the near field terms ^ for 1 < m < n and for any (n, q) G 

N" X N, 

■^^{Uiq,-)=ai_^,q,m,-, (5.14) 

and, for any (n, q, to) G such that, n> 1, and 1 < to < rr. 


L"/3J 


^-m{ui,q) = - E 


n—3k,q, — m-\-3k 


^_J2C-m+3k,2k,^A^nSy + (5.15) 


k^l 


2=0 


fe=l 
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Here, dn,q,m are defined by 


Ln/3J fc 

V(n,<7, /)GN X Z, Qn,( 3 ,m,- = ^ ' Qw-3fc,q,m+3fc,- ^ ' C*m+ 3 fc, 2 fc,i,-(In (5) , 

fc—0 i—0 

with 

n j 

V (j,m) e N X Z, _ = C(so,j), and ^ ^ XIif > 0- 

± fc=max(l, —m) 

5.3 Construction of the terms of the asymptotic expansions 

The matching conditions then allow us to construct the far field terms ^ and 11^ and the near field terms 
C/^ q by induction on n. The base case is obvious since we have seen that the macroscopic terms Uq ^ are entirely 
determined by Proposition (3.10), the boundary layer correctors Hq q are defined by (2.30), and the near field terms 
=0,geNby(4.22). 

Then, assuming that ^ and ^ are constructed for any m < n — 1, we will see that (5.11),(5.13) and 
(5.14),(5.15) permit to define both ^ (and consequently H^ g) and ^ j. for any q € N. 

Far field terms We remind that, for a given g € N, the complete definition of the macroscopic terms ^ only 
requires the knowledge of the for any integer m between 1 and n. In fact, the conditions (5.13) define 

exactly g): in the right-hand side of (5.13), the quantities known is uniquely defined) 

and ^ _i_) are known since ^ _i_ are already defined (induction hypothesis). In addition, since 

{ n—3fc <J 

X X^-^(“"-3kp)f-m-r3/c(S-r,9-p) if" S/c 7 ^ 0 , 

± r=max(l,m-3fc) P=0 

tm+3kiso,q) ifn-3/c = 0, 

the coefficient q -m+ 3 fc -i- i® "'^ii defined for any k such that 1 < fc < [n/3j (£lij{un- 3 k,q) is known by the 

induction hypothesis). Naturally the conditions (5.15) define q) is the same way. Finally, the definition of 

the boundary layer terms H^ ^ follows from (2.30). 

Near field terms Similarly, the definition of the near field terms g j. requires the specification of the quantities 
■^m{Un q±) f^t any integer m between 1 and n. The condition (5.11) exactly provides this missing information 
for g Indeed, in the right-hand side of (5.11), the computation of df^_^ q m + requires the knowledge of 

{ n—m—3k q 

X Y X^-^(“"-'"-3fc.p)C+3fc('S-r.9-p) ifn-TO-3fc7^0, 

± r=0 p—0 

fm+ 3 fe(so.,) ifn-m-3fc = 0. 

for fc between 0 and [(n — m)/3j. But, since m > 0, £^ri'‘^n-m- 3 k p) ^re well defined thanks to the induction 
hypothesis. Then, (7^ ^ + is entirely determined. In the same way, the condition (5.14) allows us to define U^g _ 
as well, replacing all occurences of i^{s^^ g_p) by m € Z and all occurences £miso,q) by (so.g). 

m € Z in the previous formulas. 

Remark 5.3. We point out that the variables n and q play a very different roles in the recursive construction of the 
terms of the asymptotic expansion. Indeed, the construction is by induction only in n. At the step n, we construct 
ui,q, n^,g and Ui g^^for any qGN. 

6 Justification of the asymptotic expansion 

To finish this paper, we shall prove the convergence of the asymptotic expansion. Our main result deals with the 
convergence of the truncated macroscopic series in a domain that excludes the two corners and the periodic thin 
layer: 

Theorem 6.1. Let Nq > 0 such that 3 Nq is an integer and let denote the set of couples {n, q) € such that 
|n -f g < Nq. Furthermore, for a given number a > 0, let 

Zla = \ {—L — a, L + a) x {—a, a). 
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Then, there exist a constant i5o > 0, a constant C = C{a, Sg) > 0, and a constant k = k^Ng) > 0 such that for 
any S G (0, (5o) 

||^^_ ^ < C'<5^‘>+5(ln<5)'=. (6.1) 

(n,g)Gr>No 

Remark 6.2. With more sophisticated techniques than applied in this article it is possible to prove that the power 
of\w5 in the previous theorem is k = k{Ng) = [^(iVo + |)J. The first logarithmic term appears for Ng = 


6.1 The overall procedure 

As usual for this kind of work (See e. g. [24] (Sect. 3), [21] (Sect. 5.1), [18] (Sect. 4)), the proof of the previous 
result is based on the construction of an approximation of the solution on the whole domain obtained 
from the four tmncated series (at order Ng) of the macroscopic terms, the boundary layer terms and the near field 
terms: 


- The truncated macroscopic series Nq ■ introduce the macroscopic cut-off function 



which is equal to 1 for I a; 11 > L, and which coincides with x(^) in the region |a:i| < L —15 (see Fig. 5). The 
cut-off functions x, x± nnd Xmacro,± are defined in (1.8)-(1.17) and (4.7). We then define the macroscopic 
approximation as follows: 

uLcvo,No = xLcro.total(x) 

{n,q)^Do 


- The truncated periodic corrector series : it is given by 

n^„(x) = (i-x(x2))x+(^^k-(^^) E (6.4) 

{n,q)^Do 


The function x+{ ^^s permits us to localize the function (x) in the domain \xi\ < L while 

the introduction of the function (1 — x{^ 2 )) ensures that (x) vanishes on Td. 

- The truncated near field series ±(x): 


<.±W= E 

{n,q)^Do 


(6.5) 


We shall construct a global approximation that coincides with ^ in the vicinity of the two corners, with 
in the vicinity of the periodic layer and with away from the corners and the periodic layer. To do so, 

we introduce the cut off functions 

= ““ 4(x) = x(^). ( 6 . 6 ) 



Figure 5: Schematic representation of the cut-off function xlacro total 
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where r]{6) : K’*' —)■ K’*' is a smooth function such that 


/ 

lim r](S) = 0 and lim —— = +oo. 

(5^0 <5^0 6 


(6.7) 


For instance for s G (0,1), r]{6) = (5® satisfies these conditions. Finally the global approximation of is 
defined by 

MWo(x) = X+(x)t/jVo,+ + X-{x)Uno,- + (l-xt(x)-X-(x))(<acro.Aro(2^) + nWo(a:))- (6.8) 

Note that belongs to but does not satisfy homogeneous Neumann boundary conditions on F^. 


The aim of this part is to estimate the iF^-norm of the error in (We remind that G 

is the ’exact’ solution, i. e. the solution of Problem (1.3)). It is in fact sufficient to estimate the residue and the 
Neumann trace c)„e^^. Then, the estimation of ||e^^||directly results from a straightforward modification 
of the uniform stability estimate (1.5) (Proposition 1.1): there exists a constant C > 0 such that, for S small 
enough, 

l|eArollHi(n'5) < 67 (||Ae^^||i2(QS) + ||c)„e^|^||i2(r6)) . (6.9) 

The main work of this part consists in proving the following proposition: 

Proposition 6.3. There exist a constant C > 0 and a constant (5o > 0 such that, for any e > 0, for any 5 G (0, 5if), 

llAe'^JU2(o.) + ||0„e^JU2(r.) < 67 V . (6.10) 

As a direct corollary, choosing ri{5) = sf6, e = we obtain the following global error estimate: there exist a 
constant C > 0 and a constant Jq > 0 such that, for any 6 G (0, Sq), 

r A^’o 5 

lleiVollsHf^'’) — 67(5 2 2. (6.11) 

Since coincides with — J2{n q in 62 q, for S small enough. Theorem 6.1 follows from (6.11) 

and the triangular inequality. 


The remainder of this section is dedicated to the proof of Proposition 6.3. Although long and technical, the proof 
is rather standard. 

Remark 6.4. We emphasize that is certainly not the best choice to minimize the global error. As shown 
in [16], a global estimate based on the truncated far and near field terms obtained by the compound method might 
provide a better global error. Nevertheless, for the sake of simplicity and since we are mainly interested in the 
macroscopic error estimate (that can always be made optimal thanks to the triangular inequality), we prefer using 
here 


6.2 Decomposition of the residue into a modeling error and a matching error 


Remarking that the supports of the derivatives of x% and xl are disjoint (for (5 small enough), using additionaly 
that A(7^^ ^ = 0 , we can see that 


-Ae^o — 6 ^mod “f 6 ^match; 


( 6 . 12 ) 


where, 

6^match = -^[A,X±(x)]((77Vo,± - ^acro.ATo ( 2 ^) “ (6.13) 

± 


and 

6’mod = /-(!- X+(x) - X-(x))A (MLcro.Aro(2;) + ^Noi^)) ' (6-14) 

Here, f match represents the matching error. Its support, which coincides with the union of the supports of Vx+ 
and Vxl, is included in the union of the rings ri{S) < |r=*=| < 2ri{S). It measures the mismatch between the far 
and near field truncated expansions in the matching zones. Smod, representing the modeling error (or consistency 
error), measures how the expansion fails to satisfies the original Laplace problem. 


Similarly, is it easily seen that 

= Yj 9nxi{^) (UnoA - n^o(^)) ’ 

± 

so that the error on the boundary data in supported in the matching areas. Therefore its treatment will be similar to 
the one of the matching error. 


In the next two sections, we shall estimate in turn £^mod (Section 6.3, Proposition 6.5) and fmatch (Section 6.4 
Proposition 6 . 8 ). The proof of Proposition 6.3 results from (6.12)-(6.15) and a direct application of these two 
propositions. 
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6.3 Estimation of the modeling error 

The present section is dedicated to the proof of the following error estimate: 

Proposition 6.5. There exist a positive constant Cmod > 0 and a positive number (5o > 0, such that, for any 
S G (0,(5o), 

/ S \^° 

\\^mod\\L 2 {nS-j < Cmod J ■ (6.16) 

We first note that the intersection of the supports of VxLcro,total (x) (and AxLcro.totai (x)) and 1 - xt (x) - xi (x) 
is included in the set _ 

= {(a;i,a; 2 ) G \xi\<L - y^rjiS) - 4S^}. 

Moreover, xLcro.totai(x) = x(x) on this set. As a result, 


/-(l-xt(x)-X-(x))AMLcro.tVo(a;) =-(l-X+(x)-X-(x))lf2^^jA,x(y)] I 

y(n,9)er>No 

(6.17) 

where 1^5 ^ denotes the indicator function of In the previous formula, we used the macroscopic equa¬ 

tions (1.13) (the functions ^ are harmonic in fix U CIb unless for n = q = 0 where — Auo,o = 0). On the other 
hand, 

-(1 - xtW - xl(x))A<„(a;) = -(1 - xt(x) - xl(x))la.^^ 5i"+«An^,, + (6.18) 

{n,q)GDNg 

where. 



{n,q)GDNg 

+ (1 - xt(x) - xl(x))[A, x+(^^)x-(^^)(l - x{x,))] E (6.19) 

(n,g)er>jVo 

Collecting (6.17) and (6.18), we end up with 


^mod — ^mod.l “f ^mod,2 7 


( 6 . 20 ) 


with 

S.od .2 = -(1 - xl(x) - Xl(x))lni, E + [A,x(y)]<,(x)) . (6.21) 

(n,q)GD]Vg 


■^^mod.! is supported in a domain where \x 2 \ > Cp{6) wherein the periodic correctors ^ are exponentially 
decaying. As a result, iflmod,! converges super-algebraically fast to 0. More precisely, we can prove the following 
lemma, whose proof is left to the reader: 


Lemma 6.6. For any A G N, there exists a positive constant Cn > 0 and a positive number So > 0, such that, 
for any 6 < Jq, 

\\£mod,lh2^a^)<CN6^. ( 6 . 22 ) 

It remains to estimate £mod, 2 - Analogously to the case of an infinite thin periodic layer, we naturally use the periodic 
correctors equations (1.15). Nevertheless, the estimation requires a careful analysis because the fields 11^ ^ and 
uf^ g are singular. We prove the following lemma, whose proof is postponed in Appendix F. 1: 


Lemma 6.7. For any £ > 0 sufficiently small, there exists a positive constant C > 0 and a positive number Jq > 0, 
such that, for any 5 < Jq, 


||^mO(af,2 II — ^7 



(6.23) 


Obviously, Proposition 6.5 is a straightforward consequence of (6.20), Lemma 6.6, and Lemma 6.7. 
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6.4 Estimation of the matching error 

We now turn to the estimation of the matching error: 

Proposition 6.8. There exist a positive constant Cmatch > 0 and a positive number Sq > 0, such that, for any 

S < (5o. 


and 


II II(H'^) — ^match 

' J \ -^o-e-5/3 




ll^n^Tyj^ ||£,2('p6^ Cffiatch I ^ 


p{6) 


+ 5 ^p{S) 


No-i-e 


(6.24) 


(6.25) 


We shall evaluate ||Cmatch||L2(f25) and ||9ne^^ llL2(r5) in turn. Both of these functions are supported in the overlap¬ 
ping areas. The function fmatch is supported in the overlapping areas fljJiatch 

^match = {{xi,X2) G ^^,v{S) < < 2p{6)} . 

and 11 On 1 1^2 (pi) is supported in 


r^ch = {{xi,X2) G T\p{S) <r^ < 2p{6)} . 


We shall 


estimate ||£^match|!p 2 (Q+^^j (resp. ||(9„e^J|p2(r+__^j)butasimilaranalysis 


can be made 


for ||£^match|li,2(f2-^^_^; 


(resp. l|a„e^J|p,(p-_^^))). 


We start with the computation of Uff^ _|_ — Aro(®) ~ which, thanks to the matching conditions (5.10), is 
expected to be small. The following computation is based on an expansion of the truncated series of far and near 
field terms in the overlapping area. 

6.4.1 Expansion of w^acro Afo ’^bl,No and + in the overlapping area 

For any couple (n, q) G we consider the integer k{n, q, Nq) given by 

I I {Nq — q) — n Nq and q have the same parity, 

k{n,q,NQ) = l . 

I %{NQ — q) — n—k otherwise. 


Macroscopic truncated series u^acrow • ^n view of (3.28), for any {n,q) G there exists a function 

7^™cro.n.,.^ G V^ 2 '^(F!t) H for any /3 > 1 - 2(fc(».9Wo)+i) 

q ^{No+r-q)-n 

^n,q ~ ^n,q '^niacro,n,q,S ^n,q ~ ^ ^ r,m,+(^~*") ^ ) In r~*" ). (6.26) 

r=0 m——n 

The reader may verify that, for non negative integer r, k{n, 9 , iVo) + LtJ — Li i^o + r — q) — n\. Then, a direct 
computation shows that 

^ d^-^+'ihijR+S,9+) = 

{n,q)GDNg 

n ^{oin,q+m) 

(n,g)e-DNo m=-an,q r=0 

where a„g = ^(NQ — q) — n. Then, using Lemma 5.1, reproducing the calculations of (5.6) and using the matching 
conditions (5.10), we see that 


§(an,q+m) 


^ <5^+«<,(i?+<5,0+) = 

,5§"+9 

E 

V 

iR+)i"^-^Wm,r,+ {0+,\nR+). 

(n,q)eD]Vg 

(n,q)eDjvg ™ 

— — dr. 

l,g r —0 

(6.27) 

Finally, noticing that xLcro,total(a^) 

- L 

Amacro,+ V (5 

X2 \ 
’ <5 / 

in truncated macroscopic series 

can be written as 

/ 


n 

1(0.1,9+™) 

\ 

^macro,Afo ~Xmacro,+ (X ) [ 

si"+9 

E 

V 

^n,q,m,+ / , 

(i?+)i™-’'u;„,,,+ (0+,lni?+) 

V 

(n,q)eDjvg 

— — dn 

,q r=0 

) 

+ '^Lcro.Aro(^) (6-28) 
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where 


'^macro,Afo “ X™acro,+(X^) ^ ~*'®’7?-macro,n,g,(5 (x) . 

(n,q)GDNg 


(6.29) 


Boundary layer correctors series Similarly, the asymptotic formula (3.31) for the periodic corrector 

g associated with the matching conditions (5.10) gives 

/ n |(a„,,+m) 

nk^„(x)=X-(X+) ^ ^ ^ |X+|i— >„,,,+ (ln|X+|,X+) 


Jn,q)eDis,g 


m— — OCn 


r =0 


+ ^k,tv„(x). (6.30) 


The remainder 7^^^. Nq (^) written as 


K 


^k.Wo(x) =X-(X+) ^S”+^7^^,acro.„.,.^(x), (x) = (xi, 0)) (X) 

{n,q)eDisig j=0 

(6.31) 

where the functions H for any /3 > 1 — ^ W„^qj G 12+(S) and K 

is a given integer depending on n and g. 


Near field truncated series The derivation of the expansion of the truncated near field expansion 

is much more direct. It may be directly obtained using formula (4.24) (taking K = [a„,qj): 

( n |(a„,,+m) 

(n,9)er>jVo m=-a„^g r=0 

( n |(a„,<,+m) 

^ ^§"+9 ^ J2 |X+|t'"->™,..+ (ln|X+|,X+) 

{n,q)GDNo m=-a„,g 


No,+ 


r—0 


+TZ' 


NF,No^ 


where 


(6.32) 

'^%F,No = 'J®”+®7^^^F.n.,(X) (6.33) 

{n,q)GDNo 

T^nF. n,qO^) belonging to 23^ .^(^ 2 "'’) for any /3 < 1 + ^ g ( 1 / 2 , 1 ), 7 — 1 sufficiently small. 


6.4.2 Evaluation of the remainder 

Subtracting (6.28) and (6.30) to (6.32) gives 


UnF,No ~ "“macro,ATo ~ If^o — '^NF,Ng ~ T^BL,No ~ ^macro.TVo' 


(6.34) 


To evaluate the matching error, we shall consider separately the three terms of the right hand side of the previous 
equality, estimating their norm and the norm of their gradient over the domain proof of the 

following three lemmas can be found in Appendix F.2. 

Lemma 6.9 (Estimation of the macroscopic matching remainder). For any e > 0, there is a positive constant 
C > 0 such that 

\\<acro,No\\LHni^j<Cp{S)^°-^+l, and (6.35) 

Lemma 6.10 (Estimation of the periodic corrector remainder). For any £ > 0, there is a positive constant C > 0 
such that 

I1'^bl,jvoIIl2(s2+^^j < 77(5)^«+3-= \\'^T^BL,No\\LAn+,j + \\T^BL,No\\LAdQ+^,^,nr^) ^ C'(5"1/277((5)^°+3-^ 

(6.36) 

Lemma 6.11 (Estimation of the near field matching remainder). For any £ > 0, there is a positive constant C > 0 
such that 


\'^NF,No\\m<:il,J <C6(^ 


p{S) 


No-e-2 


and 


ll'^AfF.7VollL2(an+,^^nr3) < 


nf,NoI\la^iLJ - ^ (^) 

^r 13 

5 


A^o-e-l 


p{6) 


(6.37) 


(6.38) 
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Finally, the proof of Proposition 6.8 is a direct consequence of (6.34), Lemma 6.9, Lemma 6 .10 and Lemma 6 .11, 
and the following inequalities: 


||Vx+(x)||^„„(f 2 + ) < (777(5) \ ||Ax+(? 

^ match/ 


< Cr]{sy 


A Far and near field equations: procedure of separation of the scales 

In this appendix, we explain how, using the separation of scales, we formally get the macroscopic equations (1.13, 
1.14) for the functions ^ and the the boundary layer equations (1.15) for the functions 11^ ^ with the right-hand 
side Gf^ g dehned in 1.16. 

A.l Derivation of the macroscopic equations 

Let us consider a point Xq € fix U Hb, and let us consider a vicinity V (xq) of Xg such that V (xg) C fix U rie- 
Then, if |a;gx| > L, the macroscopic expansion (1.7) gives ^(x) = ^(x), and inserting (1.6) in (1.3) 

gives 

^ = / (A.l) 

To obtain (1.13), we identify the different powers of 5^ +*, treating indexes {n + 3k, q) and {n,q + 2k) as dijferent 
powers of 5. 

If IxgxI < L, the macroscopic expansion (1.7) gives u^p j^ yx.) = X (x) 9 ( 2^11 inserting 

(1.6) in (1.3) gives 

^ 5^+«a(x(^)<,+<,(xi,|)) =/ (A.2) 

(n,(j)GN 2 

Then, we consider S such that for any x G 14(xg), \x 2 \ > 25; then the cut-off function x (x) equal to 1, and 
thanks to (1.10), we can neglect the terms AII^ ^{xi,x/S). Finally, we extract then the same powers of 5^+*, 
and we get equations (1.13). 


A.2 Derivation of the boundary layer equations 

We consider again the equation (A.2) for any |a;i | < 1 (this time, we put no restriction on X 2 ). Using (1.13), we 


E 

(n,9)eN2 



dx2 


u 


5 

n,q 



n 



An: 


n,9 



= 0 


(A.3) 


On the one hand, on the support of x'{x 2 lS), we see that X 2 is small; then we can use an inhnite Taylor expansion 
for both uf^ g and its 0 : 2 -derivative: 


00 p 

ui,q=^^9P2'^n,qixi,0^), ±X2 > 0, 

P^o 
oo p 

dx2ui^q = E ^9P+^Un,q{xi,0^), ±X2 > 0 

p=0 


(A.4) 

(A.5) 


We insert (A.4) and (A.5) in (A.3). Since the functions ^ depends on the fast variable X 2 = X 2 / 5 , we make the 
variable change in the Taylor expansions as well. Then, reordering with the same powers of 5^"'"'^, we get 


(n.gieN^ ± p =0 


+ E E X±iX2)^—^dP+\n,q-l-p{xi,Q^)^ 

± P^O 


Y, 5^+«Anl_,(xi,|) = 0 . (A. 6 ) 

(n,q')GN2 


On the other hand, expanding the Laplacian of each function 11^ separating the slow scale xi from the fast scale 
X = x/5, we get 


An: 


n,q 


(- 1 ) 


= 

^xi n,q 


(a:i,X) + -5,,5x,n^,,(xi,X) + 


Axh: 


s 

n,q 


(xi,X) 


We use this last relation in (A. 6 ) and we identify with the same powers of S‘^ 3 + 1 ^ then we get (1.15) with the 
right-hand side dehned in (1.16). 
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B Technical results associated with the analysis of the boundary layer 
problems: transmission conditions 

B.l Proof of Proposition 2.4 

To prove Proposition 2.4, it is sufficient to prove by induction that, any sequence (rtfc, 11^) solution to Problem (2.2) 
satisfies the following three properties: 


n,(xi,X) = £aP^(w,_p(a:i,0))rWp‘(X)+£ 5 ^ 71 ( 9 ,, tr,_p(xi,0))r%"(X), (B.l) 

p—0 p—1 

[w,(a:i,X)]r = ^ 9^^(w,_p(xi, 0))r + ^ 0))r, (B.2) 

p—1 p—1 

q q 

[9^2'Wg-i(a:i,X)]r = '^Af^dP^{uq-p{xi,0))r+ dP~^{doo^Uq-p{xi,0))r- (B.3) 

p^2 p^2 

Here, we posed = 0 for any negative integer p, and, forp > 0, € V^{B) is the unique decaying solution 


to 


-Ax%‘(X) = T;*(X) + ^[ 5 o(X)] + ^bi(X)] in B, 

d„Wp = Oonaflhoie, 

[dx,w^io,X2) = dx,w^ii,X2), X 2 GR, 


(B.4) 


where 


(B.5) 


f;(x) = 2ax,vp;_i(x) + ma_2(x) + (-i)Lp/2J (2(pp(x))<);''“) 

P-1 r„, ^VM P-1 

k^2 ^ k^2 

and the constants and are given by 

Vl= f F‘(X)I?(X)dX, K = - f F‘(X)V(X)dX. 


(B.6) 


In formula (B.5), is equal to the remainder of the euclidian division of p by 2 (/. e. is equal to 1 if p is odd 
and equal to 0 if p is even) and 5®''®" = 1 — (5®''®" is equal to 1 if p is even and equal to 0 otherwise). Moreover, 

[rj denotes the floor of a real number r. 

Similarly, WT = 0, forp < 0, and, forp > 1, WT e V^{B) is the unique decaying solution 


-AxlPp"(X) = f;(X) + ^[5o(X)] + •^[pi(X)] mB, 

dnWp = 0 on (9flhole, 

dx,w;{0,X2) = dx,w;{i,X2), X2eM, 


(B.7) 


where 


f;(X) = 29x,VP;_i(X) + 1T;_2(X) + (-1 )Lp/ 2J (2(5p(X))<5°dd) 


p 

P-1 


E(-1 

fc =2 


,Lfc/ 2 j [gfc(^)] 


p 
P-1 




2 -K - p—K • / , \ 2 

k^2 

and the constants and JVp are given by 

v;= [ T;(x)i?(x)dx, v; = - / F;(x)V(x)cix. 

JB Jb 


k fe+l5 


(B.8) 


(B.9) 


Remark B.l. The well posedness of Problem (B.7) and Problem (B.4) results from the application of Propo¬ 
sition 2.2. By construction, the right-hand sides of Problem (B.7) and Problem (B.4) satisfy the compatibility 
conditions {CxT) and {Cjy), ensuring that Wp and Wp belong to V^{B). 
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B.2 Base cases: p = 0 and p =1. 

The base case p = 0 has been done in Subsection 2.2.1. It is easily verified that, for p = 0, Problem (B.4) and 
Problem (2.14) coincide since I?g = Nq = 0. Besides, formula (2.13) (resp. formula (2.12)) exactly corresponds 
to (B.l) (resp. (B.2)). 

The case p = 1 has also be treated (Subsection 2.2.2). We remark that the profile function vanishes, and 
that Problems (B.7) and (2.21), which both define Wi, are similar. As a result, formula (2.20) (resp. the jump 
conditions (2.19) and (2.18)), exactly corresponds to (B.l) (resp. (B.2) and (B.3)). 

B.3 Inductive step 

Let <7 € N. We now assume that formulas (B.l), (B.2) and (B.3) hold for any non negative integer p such that 
p < q — 1. We prove that they are still valid for q. We shall follow the procedure described in Subsection 2.3. 

B.3.1 Computation of G'q(xi, X) 

This is by far the most technical step of the proof. We remind that Gq{xi, X), defined in (2.3), is given by 

G,(xi , X) = 2 {go (X)) , 0))r + ^ [go (X)] [uq (xi , 0)]r 

+ 2(gi(X))((9^,Mg_i(a;i,0))r + ^[pi(X)][(9^,'u,_i(a:i, 0)]r 

+ 2d^,dx,Tlq_, + 2dlnq_2 + Aq{xi,X) + Bqix^X), (B.IO) 

where, 

^,(xi,X) = ^2{gp{X)){dP^Uq_p{xi,0)}r, (B.ll) 

p=2 

BqixuX) =j2l[9p{^)] [dP.,Uq.p{xuO)]^. (B. 12 ) 

p=2 

We shall rewrite the different terms of the third line of (B.IO) using the following substitution rules: 

- We replace the normal derivatives of the macroscopic terms with their corresponding tangential derivatives 
using the following two formulas (we remind that for any p e N, Aup = 0 in a lower (and upper) vicinity 
of the interface T): for any non negative integers k, I and p, 

(B.14) 

- Foranyp < g, we substitute the jump of the traces [up(a;i, 0)]r (ortheirtangential derivatives 0)]r) 

for their explicit expression (B.2). 

- For any p < g — 1, we replace the jumps of the normal traces [dx 2 Up{xi, 0)]r (or their tangential derivatives 
(9^1 [dx 2 Up{xi, 0)]r) with their explicit formula (B.3). 

- We replace Ilg-i and Iiq -2 with their tensorialrepresentation (B.l). 

Computation of X) We divide the sum (B.l 1) into its even (p even) and odd (p odd) components, and 

we use formulas (B.13)-(B.14) to obtain 

X) = ^ Al,{X)dP, {uq.p{xu0))r + 0))r. (B.15) 

p^2 p^2 

where 

A‘(X)=2 (5p(X))(-l)Lp/2J<5;ven yl^(X)=2(pp(X))(-l)Lp/2J5°dd^ 
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Computation of Bq{xi, X) As previously, we first divide the sum (B.12) into its odd and even components Bq i 

Bg,2, 

Bqixi,X)^Bq,i{xi,X)+Bq,2ixuX), (B.17) 

where, using formulas (B.13)-(B.14), 

BqAxi,X) = ^ d^^^[Uq.k{x,, 0)]r, (B.18) 

k^2 ^ ^ 

and 

S,,2(xi,X) = ^ ('(-1)L'=/2J dl-\d^,Uq.^(x^MT- (B.19) 

We shall evaluate and Bq^^ in turn, using the induction hypotheses (B.2)-(B.3). In view of (B.2), we have, 

k—2p—1 

k—2 p—1 

We remark that the last term of each sum over k, corresponding to k = q, can be removed (since the inner sum is 
empty in this case). Then, the change of index p ^ p + k gives 

BqpixuX) = i?‘,i(X) dl{uq.p{xi,0))r + S"i(X) d^,;\d,,Uq.p{xi,0))r, (B.20) 

p—3 p—3 

with 


p-i 

i3‘,i(x) = y](-i) 

fc=2 


[k/2j 


k(x)] 


P-1 


^even 


7^1 

^p—k 


and 


i3;,i(x) = y](-i) 


L/C/2J 


bfc(x)] 


^even 


^p—k' 


(B.21) 


fc=2 


Similarly, thanks to formula (B.3), 

^,.2(X1,X) = [5fc(X)] 9P+'=(u,_fc_p(xi, 0))r + 

k—2 p—1 

k—2 p—1 

which, using the change of index p ^ p + k, yields 

Bq^xuX) = E Sp. 2 (X) dl{uq.p{xi,0))r + E 5p"2(X) 0))r, (B.22) 

p—3 p—3 

Here, 


p-i 


^p.2(X)=E(-1) 


Lfc/2J 


[fffe(X)] 




P-1 


p —fc + 1 


and i?" 2 (X) = E(-l) 


lk/2j 




fc=2 


fe=2 


;Codd \ rn 

2 ^k -^^p-fc+1- 


(B.23) 


Finally the sum of (B.20) and (B.22) leads to 


BqixuX) = EK,i(X) + i?‘,2(X))5^^(u,_p(xi,0))r 

p=2 

+ E(^p" i(X) + ^p"2(X)) d^,;\d,,Uq.p{x,,0))r, (B.24) 

p=2 

Here, we have artificially added the term corresponding to p = 2 in the two summations, using the convention that 
the constants i ?2 i and ^2 2 2 vanish (in the definitions (B.21) and (B.23), the sums are empty). 
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Computation of 2dx^dxi^q-i and dx^Ilq -2 These computations are less technical. The differentiation of 
formula (B. 1) (recursive hypothesis on the tensorial representation of Ilg) with respect to both xi and Xi gives 


q-l q-1 

= J2dPt\uq.i.p{xi,0))rdx,W^iX)+J2d^xAdx,Uq.i.p{xu0))rdx,W;iX). 

p—0 p—1 


Then, making the change of index p p 1 and using the fact that dxi Wq = 0 (see formula (2.15)), we get 


q 

2dx,dx,Il,_,{x„X) =J2dP^{u,_p{xi,0))r {2dx,W^_,m) 

p^2 

q 

+ J29Vi9.,Uq.p{xuO))r {2dx,W^.AX)) . (B.25) 

p=2 

Analogously (differentiating formula (B. 1) twice with respect to xi, then making the change of index p ^ p + 2), 
dlXq-2{xi,X) = ^5^, 0))r%‘_2(X) 5 ^ 71 ( 4 . w,_p(xi,0))r%"_2(X). (B.26) 

p=2 p=2 

Here, we have use the fact that vanishes. 


Summary Collecting formulas (B.IO), (B.15), (B.24), (B.25), (B.26), we end up with 

G,(xi,X) = 2(go(X))(u,(a;i,0))r + ^[5o(X)]K(a:i,0)]r + 2(pi(X))(a,,zt,_i)r 

+ \ [ 5 i(X)] 0)]r + ^ Fl{X)dl^ {uq.^{x^,Q))T + ^ F;{X)dlp{dx,Uq.^{x^,G))v, (B.27) 

p^2 p^2 

where, 

f;(X) = a;,(x) + i?‘,i(X) + Bl^{X) + 2ax,%'_i(X) + %‘_2(X), 

and 

f;(x) = A^(X) + BI,{X) + BI^{X) + 2dx,w;_,x) + w;_,(k). 

Here, the functions A*, AjJ are defined in (B.16), the functions i, B^ i in (B.21) and the functions B^ i, B^ 2 
in (B.23). Of course, it is easily verified that the preceding definitions of Fp and F^ coincide with the ones given 
in(B.5)and (B.8). 

It is important to note that, for any fixed xi € (—L,L), Gq(xi, •) belongs to (12“ (,8))^ because it is the combination 
of exponentially decaying terms (more precisely, functions or first derivative of functions belonging to V^(B) and 
compactly supported ones). 


B.3.2 Computation of the normal jump [dx 2 Uq-i{xi,Q)]T 

As in Subsections 2.2.1, 2.2.2 and 2.2.3 (base steps), for any fixed xi € {—L,L), Proposition 2.2 ensures the 
existence of a periodic corrector I{q{xi, •) € )2+(8) satisfying (2.2-right) if and only if Gq{xi, •) is orthogonal to 
both N and T) (that is to say Gq satisfies the two compatibility conditions (Cx))-(C 7 >/)). In view of the second and 
fourth lines of Lemma 2.3, we have 


9 r 

Gq{xuX)M{X)dX=[dx,UqXv + y^dl^{uq_p{xuQ))r / F^{X)MiX)dX, 

' p=2 

q f 

Fy,dlp{dx,Uq.p{x^,Q))T / F;{X)M{X)dX, (B.28) 

JB 


In the previous formula, we recognize the constants A/jJ and defined in (B.6)-(B.9). Finally, the compatibility 
condition {Cj\f) is fulfilled if and only if 

<? 9 

[dx:,Uq-l]T = ^dP^{Uq-p{xi,Q))TMp {dx^Uq-p^Xl,^))^^. (B.29) 

p=2 p=2 

and formula (B.3) is proved. 
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B.3.3 Computation of the jump , 0)]r 

Similarly, the first and third lines of Lemma 2.3 give 


[ G,(xi,X)I?(X)dX = -[M,(xi,0)]r + [ 2{g,{X))V{X)dX+ 

JB JB 

^aP^K_p(a:i,0))r [ F^iX)ViX)dX, + ^dP;^{d,,u,.p{x,,0))r [ F;(X)I?(X)dX, 

..,-9 JB „_o JB 


Here again, we recognize the constants 22* and defined in (B.6)-(B.9). If we additionaly remind that i^"(X) = 
2{gi{X)) (cf. (B. 8 )) and that VL* = 0 (see B.5), the compatibility condition (Cp) is fulfilled if and only if 


K(a;i,0)]r = ^X>‘aP^(Mq_p(a;i,0))r+ ^22^9^^ ^{d^^Uq-p{xi,Qi))T 

p—1 p—1 

which ends the proof of formula (B.2). 


B.3.4 Tensorial representation of H, 

Under the conditions (B.3) and (B.2), Problem (2.2-right) has a unique solution in Moreover, substituting 

[Mq(a;i, 0)]r and [(9a;jMq_i(a;i, 0)]r for (B.2) and (B.3) in formula (B.27), using the linearity of Problem (2.2-right), 
and reminding the definitions of Wp and Wp (B.4)-(B.7), we see that 

n,(a:i,X) = ^aP^(u,_p(a:i,0))rlUp*(X)+ ^aP-*(a,,w,_p(a:i,0))rVL;(X). 

p—0 p—1 

which is nothing but formula (B.l). 


C Technical results associated with Analysis of the macroscopic problems 
(macroscopic singularities): proof of Proposition 3.8 

The proof of the Proposition 3.8 is by induction. Before starting the proof, we first need to present a technical 
Lemma which will be usefull to define the functions Wn,p,±- 


C.l A prelimininary technical Lemma 

We consider the infinite cone K, of angle 3tt/ 2 (i. e. in 2D, K, is an infinite angular sector of angle 37 r/ 2 ). 

/C = /(rcos0,rsin0) € r > 0, 0 € (0, ) (C.l) 


that we divide into two disjoint cones /C* and of openings 7 and ^ — 7 , 0 < 7 < 

/C* = {(rcos0,rsin6l) G r > 0, 0 G I*} = {(rcos0,rsin6l) G r > 0, 0 € 


(C.2) 


where I* = (0, 7 ) and = ( 7 , ^). We denote by L.^ the interface between /C* and /C^, i. e., L.^ = O/C* n 9/C^. 
In the cases we are interested in the cone K, is given by /C + or /C - (defined by (3.4)), so that 7 = f (for /C“) or 
7 = TT (for IC~^). In the upcoming proof, we shall use the following technical Lemma, which provides an explicit 
formula for a transmission problem in a cone for a particular right-hand side (we remind that A defined in (3.11) 
denotes the set of singular exponents). 

Lemma C.l. Let g G N, A G M, (a, b) G and 


N = 


q 

<^4-1 if X A. 


There exist N -\- 1 functions gp G C°°{X^) H C°°iX'^), (0 < p < N), such that the function 

Av = 0 in K} n K?, 


r(r, 0) = 



lnr^(/p( 0 ) I satisfies 


v(r, 0 ) = 0 , 

0 , 


Stt 

v(r,-) = 


(C.3) 


[vjr., = ar 


lnr« 


[i9ev]r = br^ lnr«. 
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We refer the reader to the Chapter 3 in [28] and the Section 6.4.2 in [26] for the proof. 

Remark C.2. It is possible to show that the functions Pp are uniquely determined if X ^ A. If X S A, the functions 
Qp are uniquely determined except for p = 0. However, the uniqueness may be restored by imposing the following 
orthogonality condition: 

/ 9o{d)wx{0) = 0 wx{0) = sm{Xe). 

Jo 

Remark C.3. Note that the set of singular exponents X (for which we have to take N = q + 1 j coincide with 
the singular exponents of the problem without transmission condition on r..j,. In particular, these exponents do not 
depend on the location of the interface F-,,. 

C.2 Recursive definition of the families Wn,p,± 

Based on Lemma C.l, we shall be able to define the set of functions G G N} of the Proposi¬ 

tion 3.8. We begin with the construction of Wn,p,+- We remind that the functions Wn,p.+ are polynomial with 
respect to In r+, i. e.. 


Stt, 


W; 


i,p,+ (6i+,lnr+) = ^Wn,p,q,+ {0^){\vir^y, n € If ,p G N, 'Wn,p,q,+ G C~([0,tt]) nC°°([7r, —]) 


9=0 


Their construction is by induction on p. The function o,-i- has already been defined in Proposition 3.2: 

Wn.o,-i-(6''‘',lnr+) = sin(A„6»+). 

Forp > 1, we construct of the form (C.4), such that the function 

v„,p^+(r+,6>+) = (r+)^"“2’w„,p,+ (6'+,lnr+) 


(C.4) 


(C.5) 


satishes 


Av„,p,+ = 0in/C+’^ n/C+’^ 

/Stt _ 

''^n,p,+ { 2 ( ~ 

[v„,p,+ (r+,7r)]aK+.inaK:+.2 = (r+)^""P a„,p,+ (lnr+), 

^ [9e+v„,p,+ (r+,7r)]aK:+.inaK:+.2 = (r+)^""Pb„,p,+ (lnr+), 


VnGZ*, VpeN* 


(C.6) 


where 


= {(r'*" cos 9'^, r'^ sin^"*") e /C*", 0 <0+ < tt} cos0^, sin0+) S /C*", tt < 0''‘ < 


and 


p-i 


^71,p,- 


^ ^ (^p—r5n,r,j)—r,+^ '^p—r9n,r,p—r,-\-(}^^ )) ? 

r—0 

p-1 

b„,p.+ (lnr+) = + A/'p'Vi_>" ,,_p_^ ^(lnr+)) . 


(C.7) 

(C.8) 

(C.9) 


r—0 


The functions gn^r,q,+’ 9n,r,q,+’ hh^r,q,+’ ^n,r,q,+ defined by the following relations: for n S Z*, r G N, g G N, 


d<! 


(p+)A„ r 9 [(p+)>" ^Wn,r,+ {Tr,lnr+))g,c+-^ndK+-A 


(C.IO) 


gq 


-1 


(p+)A„ r 9 ^ ('_^^9_____ [(p+)A„ r 1 (0^+(tt, Inr+))aK:+.inaK:+.2] , (? > 1) 

(C.ll) 

('p+)A„-r-9 hYg^^i\nr+) = ^^(~l)'^ (gp+) 9 +l {Wu,r,+ {T^Mr^))dlC+Jnd>C+.^] 


gq 


(p+)A„ r q hy^^^ j^(\nr+) = r+{-lY [(r+)^" ’’ ^ {de+Wn,rA'^M'r""))dK+JndK+A 


(C.12) 

(C.13) 


and fln,^,o,+ = 0. 
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The recursive procedure to construct the terms tUn,p,+ is the following: assume that the terms Wn,q,+ are known 
for any n G Z* and g < p — 1. To construct Wn,p,+ we need to know the source terms a„ p + and b„ p+, which 
require the computation of 5 "^, p_r + for r < p — 1. But, p* p, p_r ^ and p"_r,p-r,+ only depend on 

the function Wn,r,+, which, thanks to the induction hypothesis is known. Similarly, to construct b„ p+, we need 
to define ^ and /i" ^ ^ for r < p — 1. These terms only depend on Wn^r,+, and, as a consequence, are 

known. Lemma C.l then ensures the existence of Wn,p,+- 

Of course, we could have written a recursive formula to obtain g^r q + (resp. gnr q +’^ri r q r 9 +) starting 
form r 0 + (rcsp. p" ^ 0 r 0 r 0 +) shall not need to explicit this formula. However, we can 

prove the following useful relation: 

^n,r, 9 ,+ 9n,r^q+l,+ and ^n,r, 9 ,+ 5n,r,9+l,+ ’ (C.14) 

Remark C.4. IfXn — p € A (which is the case as soon as p is even, except if Xn — p = 0), the functions Wn,p,+ is 
not uniquely defined by (C. 6 ). In that case, we add the somehow arbitrary condition 


w^n,p.o,+ (^'"^.lnr+)w„,o.+ (6'’^)rf6»+ = 0 


(C.15) 


to restore the uniqueness (see (C.A) for the definition ofWn,p,o,+)- As a consequence, we can see that the sum over 
q in (C.4) goes from 0 to [p/2j (and not p): in other words, Wn,p,+ is a polynomial of degree at most [p/2j. 


Similarly, we can define by induction the functions Wn,p,- of the form 

b/2j 

Wn,p,-iO~ ,\nr~) = ^ Wn,p,q-iO~){\nr~y nGZ*,peN, u;„,p,g _ G 0]) n C°°([0, tt]). 

9=0 

(C.16) 

Again, the function Wn,o- has already been defined in Proposition 3.2: 


Wn,o,-{0 ,lnr )=sm(A„(6» - 7r/2)), 

Forp> 1, we construct Wn,p- such that v„^p__ = Inr“) satisfies 

lS.Vn,p- = 0 inn 


(C.17) 


'^'ri,p,-{ 2 ^ ~ 


< v„,p,_( 7 r) = 0, VnGZ*, VpGN*. (C.18) 

[•Vn,p,-ir~,0)]oK-xndK--^ = a„,p _(lnr"), 

[ [9e-v„,p_(r",0)]aK:-.inOK:-.2 = b„,p _(lnr”), 


Here, 


/C 4 _ I ("p (jQg Q ^ p gj]^ P ) G /C , —TT < 0 <0}/C ’^ = {{r cos 9 ,r sin0 )g/C ,O<0 <7r}, 

(C.19) 

and 


p-i 


a^.p,—(lar ) ^ ^ (^p—rffn.r.p—r,—^ ) A 72p_7.Pyjp. p—,. _ (lu T 

r—0 
p-i 

hn,p,—(lar ) ^ ^ {f^p-\-i—r^n,r,p—r,— i^^t' ) T p_y. _ (In r 


(C.20) 

(C.21) 


r—0 


The functions gn r q 9n r q Ki r q - ^ ^nrq- ^re defined by the following relations for n G Z \ {0}, r G N, 
p G N, 

B‘t 

p*,^,,,_(lnr-) = {wn,r,-iOMr-))a!c+xndK+.^] 


-1 


(p-)A„-r-9 gl^^^^ _(\nr-) = {de-Wn,r,-{0, \nr-))9,c-.^ndK--A {q > 1) 


(r )^" )=' 


59+1 
(9r“)9+i 
di 


(p )A„ ))a/c-.in3K:-.G 


(r )^" ^".r. 9 .-(lnr )=r [i^ )^" ^ (ae-'u;„,^_(0,lnr ))9ic-xndK-.^] 
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and gn ro - = 0 - we have; 


'^n,r,q,— yn,r,g+l,— '^n,r,q,— yn,r,q-\-l.—' 


(C.22) 


To restore the uniqueness of Wn,p,- if A„ — p € A (i. e. p even), we impose the orthogonality condition (see (C.16) 
for the definition of Wn,p,o,-y) 


/ TT 

'Wn,p,O-{d~)Wn,O-{0~)de~ = 0 . 


(C.23) 


C.3 Construction of the family g of Proposition 3.8 

We are now in a position to start the heart of the induction proof of Proposition 3.8. The proof of the base case 
and the proof of the inductive case essentially rely on the same arguments, the inductive step being as usual more 
technical. 


C.3.1 Base case 

I satisfies the following problem 

_As+ — 

o+ = 

* —m,l 

[sl„,i(a:i,0)]r = 

, [a^2sl„,i(a:i,0)]r = 


0 in n 
0 on r^), 


'D'i (sl„ o(a:i, 0 ))r + T)” {d^^st^^o{xi,0))r, 

^2 0 ))r + A^" 0(2^1, 0 ))r- 


(C.24) 


We first observe that the right-hand side of (C.24) is singular (see the asymptotic (C.25)-(C.26) below). As a 
consequence. Proposition 3.1 does not apply, at least directly. However, we can write an explicit asymptotic repre¬ 
sentation of the jump values of s+ ^ in the vicinity of the two comers. Then, we shall see that we are able to make 
an appropriate lift of the singular part of the jumps to reduce the problem to a variational one, which can of course 
be solved using Proposition 3.1 (or Proposition 3.2). 


Using the asymptotic expansion (3.18) of q, in the vicinity of xj, and thanks to the two formulas dx^v = 
—dr+v, 8 x 2 = valid on T, we can write an explicit expansion of 0 )]r in the vicinity of 

xj: for any fc € N*, for r+ sufficiently small 


k+l 


^ w„,o.-t( 7 r) - X> 5 'ae+w„,o.-t( 7 r)) + ri.fc,+ (r+) 


fe+1 


= 51 a„,i,+ f+(s;„_o) (r+)^"-i +ri.fc.+ (r+), (C.25) 


n— — m 


where, ri^k,+ € for any j3 > I — We used the convention ^ ^ 

for the integers n such that —m -f 1 < n < 0. a„ 1 4 . is defined in (C. 8 ). 

1/2 

Based on the same technical ingredients, we can prove that there exists a function r 2 ,k,+ in U 2 ^ (T) for any 
/3 > 1 - such that 

^5e+st„4(r+,7r) 


= E 


^n,l,+ 


(r +)^"-2 + r 2 ,fc.+ (r+), 


(C.26) 


b „4 4 being defined in (C. 8 ). Of course, a similar analysis provides an explicit expression for the jump values of 
slirn 1 close to the left corner: 


k+l 


s-m,iir 7 0)]r = E^"4-C(stm.o) ^ + r+k-ir ), 


—90-st^Ar ,0) 


where, for any j3 > 1— ^ fi.fc,- € and r 2 ^k,+V 2 'p {J^)- The constants a„ 4_4 and b„ 4 ^_ are defined 

in(C.20)and(C.21). 


fc+1 


= ^\ln,l,- ^ + ?'2.fc.-(r+), 

r n—l 


(C.27) 

(C.28) 


/3/2/ 
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Now, let us consider the problem satisfied by a function 

fc+i fc+i 

S-m.l.fc = S-m.l - (C.29) 

n— — m n—1 

defined on ffx H ffe- In this definition, we have lifted the most singular part of the jumps of ^ across F (up 
to a given order k) using the functions v„ i ± defined in (C.6)-(C.18). The cut-off functions Xl defined in the 
sentence following the definition (3.5). satisfies 

= 

5 + = 

< 

[■5-m.i,fc(a^i>0)]r = 

A direct computation shows that fk € Fl (/ is compactly supported away from the two 

corners), gk & and hk € for any P > 1 — ^ Thus, for k > 1, and choosing /? — 1 ^ A, 

Proposition (3.2) guarantees that Problem (C.24) has a unique solution belonging to iTQp^(nT) H 
that satisfies the following asymptotic close to the two corners: there exists 
P > 1 — P — 1 ^ A, such that, for sufficiently small, 

k 

straXk = ^^.31) 

n—1 

The existence of ^ is then a direct consequence of formula (C.29). Finally, the asymptotic formula for ^ 
directly follows from (C.31) and the asymptotic formulas for v„ i ± (See (C.6)-(C.18)). The uniqueness of 
follows form Remark 3.7. 

C.3.2 Inductive step 

We assume that the function p are constructed up top = g — 1 and that the asymptotic formulas (3.23)-(3.24) 
hold for p < q — 1. We shall prove that we can construct ^ G V 2 H V 2 p{Xl^), /3 > 1 + ^ + g, 

satisfying (3.22) and admitting the asymptotic expansions (3.23)-(3.24) in the vicinity of the two corners. The 
proof is almost entirely similar to the one of the base case. First, using the asymptotic formula (3.23) for q-p 
and the definitions (C.8) of a„_p+, we obtain the following formula for the jump value of q across F: for r+ 
sufficiently small. 


fk in rix n Ob, 
0 on Fb) , 

9k 1 
hk ■ 


(C.30) 


[s-m.Jr=y] y] C(sl™.,-p)a„.p.+ (r+)^"-P + ri.fe,+ (r+) 

P=1 -■m<n<{k+l)+l 

where ri ^ + G for any P > 1— . Here, we have used the convention that forp > 0, (s+ p) = 0 

if n > 0. Similarly, using the definition (C.9) of b„ p+, we see that there exists r 2 ,k,+ G for any 

P > 1 — such that, for r+ small enough 

=J2 E + ^2.fc.+ (r+). 

P—1 —m<n<(fc+l) + ^ 

Analogously, we obtain asymptotic representations of the jump values of ^ in the vicinity of the left corner: 
there exist ri_k,- belonging to V 2 ^^^(F) for any /? > 1 — and r 2 x,- in f®*" ^^^y /? > 1 ~ ^ 

such that, for r~ small enough 

9 

[■^-m.olr = ^ ' £n(s-m,q-p)3.n,p,-(f ) " ^ ~k’^l,k,- 

P=1 l<n<(fc-|-l)-|-§ 

= E E ^;!:(s-m.9-p)b„,p-(r-)^"-P-i+r 2 ,fe-. 
p=l l<n<(fe-|-l)+| 
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^n,q- and bn,g,- are defined by (C.20) and (C.21). Then, we consider the problem satisfied by the function 

q 

^—m,q,k ^—m,q ^ ^ ^ ^ p)XL ) 

P—1 —m<n<(/i:+l) + | 

Q 

Q{stm,q-p)xl^n,p,-{r~ (C-32) 

P=1 l<ra<(fc+l)+§ 

defined on fix H XIb, where we have lifted (up to order k) the most singular part of the jumps values of The 

functions v„ p ± are defined in (C. 6 )-(C. 18). The remainder of the proof is entirely similar to the one of the base 
case. 


D Technical results associated with the analysis of the near field problems 
(near field singularities) 


D.l Proof of Proposition 4.2 

Before proving the Proposition 4.2, let us first give a preliminary technical result (whose proof essentially follows 

from the fact that |fI ^ |sinf| and, for f € (— 7 r/ 2 , 7 r/ 2 ), |sinf| > 

Lemma D.l. For any € (|, ^) (/. e. < Oj, the following estimate holds: 

1 + 1 X 2+1 <p(i?+,0+) < (^ + 1 ) (l + |2f2+|). (D.l) 

We can now start the proof of Proposition 4.2 Let 7 e (^, l), p G {1, 2}, A G M and g G N. 

• Let vi = x^^^(X^)x_(X+)(i?+)^(ln(i?+))* and let us determine on which condition on /3, vi belongs 
to ..y(n+). We first note that compactly supported in the bands 1 ^ |X^| ^ 2. On this 

domain, we can bound x^^'^ by its L°°-norm. Then, in view of Lemma D.l (note that due to the cut-off 
function X-(X+), we only consider X+ < 0), p defined in (4.4) is bounded from below and from above. 
As a consequence, we have to estimate 

f f (i?+)2^(ln(i?+))29(i?+)2/3-27-2^7^+^X-+_ (D 2 ) 

4-00 Jl 


In addition, since is bounded from above and from below, i?+ is equivalent to |X+|. Moreover, 
(In4?+)^“? can be bounded by |X+|^, for any £ > 0; then estimating the integral (D.2) is equivalent to 
the estimate 

r-i 


I ^ + 12 A+2/3—2'y—2+e 


dX+ 


(D.3) 


and integral (D.3) is bounded if and only if 


2A + 2/3- 27 - 2 + e < -l,Ve > 0, 

Finally vi belongs to -^(^2+) as soon as/ 3<7 — A+|. 

• We now consider the function V 2 = x(-R+)(f?+)^(lni?+)^. Since the derivatives of x(-R+) are compactly 
supported, the only part to estimate is 

[ ^(i?+)(i?+)2/3-27-25p,0p27-4-r2p-r25^,0 |((i?+)^(Ini?+)«)}^ dX+fiX+, (D.4) 

4n+ 

for 0 ^ p ^ 2. First, we can exhibit a function (j)p{R~^) such that 

yp ((i?+)^(lni?+)«) = {R+)^-P(j)p{R+), with (/)p(i?+) = 0{{R+y), Ve > 0, 
so that we have to study the convergence of 


Jq+ 

Note that estimation of (D.5) for p = 0 and p = 1 is the same. Moreover, since the origin (0, 0) does not 
belong to the support of x, we only have to consider the behviour for large i?+. To estimate this integral, we 
split the domain into three parts: 
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- the domain Oi = $7+ n {O ^ ^ < 0} (the intersection of f2+ (which has holes) with 

a fixed width band located on both sides of the interface {{X ^< OandJf^ = 
0}): since fli is included in the band |0 ^ {X^l ^ 2, X^ < O}, the integral (D.5) restricted to fli is 
smaller than the same integral in the whole band. Then, again, we can use that p is bounded in this part 
(from above and from below), so that the integral converges if and only if 

2A + 2/3 — 2^ — 2p — 2(5p^o + —1,£>0 /3<^ — 2 P ^p,o — 

and, since 7 > 5 and p + Spp ^ 1, we can bound /3 by 1 — A. 

- the domain 172 = {2 < 1^2 \ < < 0} ('^"'o angular domains). We first note that R'^ is 

equivalent to |X^| in 172 . Indeed, + |X^|) < 2|Arj^|. Then, thanks to lemma D.l, 

p is equivalent to 1 + To evaluate the integral (D.5) restricted to 172 , we make first an explicit 
integration with respect to X 2 (note that 27 — 4 + 2p + 25pp ^ —1), then we use the equivalence of 
i?+ and \X^\ to obtain 


[ I ' (^+)2A+2/3-27-25p,o-2p+£^27-4+2p+25p,o^j^ + ^j^ + 

j —00 J 1 

<c f |7(:+|2A+2/3-27-25p,o-2p+e |27-3+2p+25p.o . (D. 6 ) 

J —00 

Since 27 — 3 + 2p + 26pp > 0, the most singular term in the right-hand side of the previous equality 
is As a consequence, 

J —00 J 1 J —00 

This integral converges as soon as/3<l — A — £/2. Finally, 


f ^j^+-^2\+2P 

J O 2 


— 2j — 26p^o — 2p+s 


-^^dX^dX^ < 


I 


converges as soon as ^ < — 1 + A. 

- the domain I 73 = 17+ n {(i?+ cos 0+, i?+ sin 0+) € | 6 *+ — 7 r| > ^}: in this domain, p is equivalent 

to i?+. Using the polar coordinates, we see that we have to estimate 



^ 2 A+ 2 ^ — 27 — 2 p—2<5p^o+e+27—4+2j?+2(5p^o dR~^ d0~^ 


(D. 8 ) 


and this integral converges if and only if 


2A + 2/3 + e-3 <-1,£ > 0 ^ /3 < 1 - A. 


This ends the proof of the second point. 

• Let u; = ■u;(X+, AT^) be a 1 periodic function with respect to 7^+ such that <+ 00 , and 

let r !3 = x_(Ar+)|Ar+|'’'“^(ln(|X+|))9r(;(A'+,X^). We have to study the convergence of the integral 

13 = [ {R+)^^-^'^-^p^^+^ (D.9) 

Jq+ 

In order to use the periodicity of w, we split the domain 17+ into an infinite sequence of similar non over¬ 
lapping domains Bn = 17+ n { —(n -f 1) < 7f+ < —n} (Bn consists of the shifted periodicity cell B, see 
(1.9)). Then, the previous integral becomes 



(^+)2/3-27-2^27+2g-|X+||;f+|2(A-l)(ij^|^+|)29 (^2(^ + 



dX. 


Note first that the inequalities (D.l) are valid in B-n (p is equivalent to (1 + |X^|) for X+ < 0). As a result 
p^'^ is equivalent to (1 -f |X4|)^+. Besides, there exist two constants (Ci, C 2 ) & (M+)^ such that 


I v +|2 

Cm <R+ < C2n{i+^-Am 


< C2n(l + |X2+|). 


(D.IO) 
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One the one hand, if 2/? — 27 — 2 > 0, using the second inequality in (D.IO) together with the fact that for 
any a G M (1 + I is bounded, we have. 




^2/3-27+2A-4-e^ 

new 


which converges as soon as 

^<^+ 7 -A-£. (D.ll) 

One the other hand, if 2/3 — 27 — 2 > 0, we can use the hrst inequality in (D.IO) and the fact that for 
any a G K, (1 + I is bounded to show that X 3 also converges under the condition (D. 11). To 

summarize, for a given parameter 7 G (A, 1), U 3 belongs to )2^ ^(0+) for any/3 G Msuchthat/3 < I+ 7 —A. 

• For the last equality, it is sufficient to investigate the two following kinds of integrals: 


14 , 1 = [ (i?+)2/5-27-2^.o^27-4+2p+25p.o^_^^+) |^+|2(A-l) (D.12) 

Jn+ 


and 

2^4,2 = 


lQ+nsupp{x'-{Xf)) 


(^+)2/3-27-25p,o^27-4+2p+25p.o |^+|2(A-1) (In |X+| )29 (X+, X2+)dX, 


(D.13) 

for a generic function w G L[q(,( 0+), 1-periodic with respect to such that ||el ^2 < + 00 . 

Indeed, by a direct computation of the partial derivatives of V 4 (up to order 2), we can see that for any 
a = (cti, 02) € such that ai + a 2 < 2, there exist three functions U4 0, U4 1 and V4 2 such that 


d^V4 = d^\d^\v4 = X-iXt)v4,2 + x'-{Xt)V4,l + X-{Xt)v4,0, 


where, for £ G {0,1, 2}, U 47 = In with A^ < A — 1, < g, and W 47 G is 

1 -periodic with respect to X^ and satisfies ||el ^2 l/^t(; 47 ||i 2 (g) < + 00 . ForZ 4 1 , repeating the arguments 
of the proof of the third point, we see that 


^4,1 < y]n2/5-27+2A-2-2i„o+e_ 

nGN 


Consequently, I 4 1 converges for any p G {0,1,2} if /3 < 7 + ^ — A — e/2. As for X 4 2 , we note that on 
the support of |Xj^| is bounded from above and from below so that is equivalent to |X^|. It 

follows that Z 4 2 converges for any /3 G M. 


D.2 Proof of Lemma 4.3 

Um,p,+ is a sum of a macroscopic contribution 

p 

Wm.p.macro = ^ (i?+) (In i?+, 0 +), 

r—0 

modulated by the cut-off function Xmacro,-i-, and a boundary layer contribution (exponentially decaying with respect 

to 1^2+I) 

P 

Km,p,BL = ^(|X4+|)^™->™,,,+ (ln(|X+|),X+,X2+). 

r—0 

modulated by the cut-off function x_ (A//"). Consequently, 

MAm,p,+ = [A,x(i?’'')](x macro, macro +X -{^1 Wm,p,BB)+X{IP')^ (X macro,BL + X-(^i+)^m.p ,macro) • 

(D.14) 

The hrst term in the right-hand side is compactly supported (since Vx(f?'*') and Ax(f?~'’) are compactly sup¬ 
ported). Therefore it belongs to .y () for any real numbers /3 and 7 . It remains to estimate the terms of the sec¬ 
ond line. The proof is technical but the main idea to hgure out is that A (xmacro,-i-^m,p,BL) and A (x- (A/" )l£m,p,macro) 
counterbalance (up to a given order). 

We start with the explicit computation of A (xmacro,-i-^m,p,macro)- We consider the function 
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already defined in (C.6). Vm,r,+ is defined in the union of the two cones and (see (C.7)). It is smooth 
on and and it satisfies Avm,r,+ = 0 in and Using the fact that Xmacro,+ = xi^}) ™ 

support of X- (^1^) (Xmacro,+ = x (^2 ) < “!)’ ^6 have 


P P 

^ (Xmacro,+ macro) = X— (^1^ ) ^ )'^ni,r, + ) + (1 ~ X— (^l")) ^ ' { ^(Xmacro,+ j )'V'm,r, + ) } ■ 

r—0 r—0 

(D.15) 

The second term in (D. 15) belongs to for any real numbers j3 and 7 because (1 — x- Vxmacro,+ 

and is (1 - X-(^i^))AXmacro, + are compactly supported and Av^ = 0. Next, 

^ {x{X^)^7n,r, + ) = ‘^x!(.X^)dx+^m,r,+ + x") Vm,r,+• 

Using the Taylor expansion with integral remainder, 


K-l 


v„,,,+ (A+,A2+) = y] 0±)^^ + 


1 




k=0 


k\ 


(K-iy. 


dtK 


{Xt -1) 


K-l 


dt, 


reminding that r,+ | is smaller than In on the support of the derivative of x(-^^)> 

it is verified that 


A(x(X2+)v„.,.+ ) = 29^+(v„,,, + ) (fffc) + [v™,,. + ] [gfc] 

fc=0 ^ 

+ (2a^-/(5x+v™,.,+)(gfc) + [gk]) + , X+), (D.16) 

I —n V ^ ^ / 


where rm,q,p is a smooth function whose support is included in the band 1 < |A^| < 2 that satisfies 

\rm,r,p\ < In |A+r, (D.17) 

for a constant C depending on m, r andp. We remind that (vm,r,+) and [v^ denote the mean and jump values 
of Vm,r,+ across the interface dJC^^ n and that the functions {gk) and [gk] are defined in (2.4). These 

functions only depend on and their support coincides with the support of xAX ^)■ To obtain (D.16), we have 
used the fact that Avm,r,+ = 0 on the support of {gk) and [gk] so that we can use the formulas (B.13)-(B.14). 
Noting that 9^+0^) = (—1)^9^+ and using the definition (C. 10) 

of 5m.r.fc,+’ we see that 

9^+(v^.,,+) = (-1)'=9^+ {((i?+)^'"-'-u;™,.,+ (lni?+))} = |A+|). (D.18) 

Analogously, since 9_y+Vm,,_|_(Aj^, 0^) = —(i?+)'’'™“’’“^9e using the definition (C.ll)of 

we have 

9^-/(9^+v™,,.+) = (DT9) 

Next, we replace the [vm,r,+] = [{R~^)^’^~^Wm,r,+] with its explicit expression (C.8), taking into account that 
Vq = Vq = 0, we obtain 

9^+K,,, + ] = (|A+|)^™-'=-'- ^VAArn.jX+r-lAAXt\)+V'^.,gAlX+r-j^^^^^ • (0-20) 

Similarly, substituting [9j^+Vm,r,+] = —(i?~'')^’"“’'“^[9e t(;m,r.+ (lni?'*',7r)] for its explicit expression (C.9), us¬ 


ing the relation (C.14) and the fact that and Afi vanish, we get 


9^-/[9;,+v^,,.+] = {\X+\A-’^-^ ^A+i-j9A,k+r-jAAXt\) +K\i-j9A,k+r-,,AAXt\)j . 

(D.21) 

Substituting the left-hand sides of the last four equalities for their right-hand sides in (D.16), summing the contri¬ 
bution (D.16) from r = 0 to p, we end up with 

^ (Xmacro,+macro) — 7 ?m,p,macro(A^l : X2 ) 
p r 

+ X-{X+) y] \X+\)Ai{XA + gAr-k,kA^^\Xt\)Al{Xt)) (D.22) 


r=0 


fc=0 
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where, 


i?m.p.macro = X-^2^ ) + (^ " X-i ’ ^2+)^™.’--+) } ’ 


r=0 


r=0 


Al{Xt) = < 5fc > + E + <5r(-l) 

r=o 


-1)L^/2J M^i 


fc—£+1 I 5 


and 


^^(X2+) = 25r(-l)L'=/2J < fffc > +E + <5r(-l)L'/^J ^A4"_,+i) • (D.25) 


(D.24) 


Comparing the previous two formulas with (B.5)-(B.8), it is recognized that 

A\ =-XWl-2dx+Wl_^-Wl_2 and Al = -- 2d- Wk_^. 


(D.26) 


Applying the first point of Proposition 4.2, we note that Rm,p,macro 
and therefore for any (3 < 2 — (Am — p). 


belongs to .^(f2+) for/3 < 7-(Am-p) + |, 


The computation of A (x- (A'('' )l4m^p, bl) (last term to evaluate in (D. 14)) is much easier. We remind that pm,r,+, 
defined in (3.32), has the following expression: 

/ r r \ 

Pm...+ (ln|A+|,X+) = E 9l,r-k,k,+ MXt\) Wt (X+) +E 9kr-k,k,+ MXt\) (X+) . 


^fe=0 


fc=l 


Then, using the definition (C. 11) of ^ ^ ^ we see that 


A{|X+|^™-’'p^,,_fc^,,+ (ln|X+|) W^,*(X+)} = |X+|^'"-'-gl,,_fc,fc^+(ln|X+|) AM/t(X+) 


+ (2pl,,_fc,,+i^+(ln|X+|)a^+VP‘(X+))+ |X+|^--’-2 (gl_,_,_fc+2,+ (ln|X+|)lC*(X+)) . 

(D.27) 


The reader may verify that a similar formula holds for A j, ^(In iXj^Dl, replacing 

i with n . It follows that 


the subindex 


A {x-{X+)U^,p,bl) = X-(X+) E + Wl 

\r=0 \fc=0 / / 

+X-(X+) (EIX+I"™-'- (E5™..-fcX.+ (ln|X+|) [AW^ + 2d^^W^_,+W^_,u]+R^,,MXt,X+), 
\r=0 \fc=0 / / 

(D.28) 

where, 

i?m,p.BL(X+,X+) = 

P+1 

X-("^1^)1 "^1^1 ^ (‘^ 9m,v-\-l-k,k.-\-(^^\^t\) ^xf^k-1 + ^5m,p+l-fc,fc,+ (^^ 

k^O 

p+2 r 

+ X-(X+) E E ('7+.-fe.fc.+(ln|++l)l^+2 +<.-fc.fc.+ (ln++|)l^+2) 

r—p+1 k—0 

+ 2Vx-(X+) • VZ7m,p.BL +Z^m.p.BLA(x-(X+)). (D.29) 

The support of the terms of the fourth line are included in a band B = supp(x(_(X(^)) n 12'*'. Since, the 
function and the derivative of 27m,p,macro are exponentially decaying in this band, it is verified that these terms 
also belong to QJ^..y(r2+) for any real numbers /3 and 7. The terms of the second and third lines are of the 
form |X(''|''’"“-P“^wi(ln |X(''|,X('',X^) or |X(''|'''™“^’“^W 2 (ln iXj^j, X+X^) where wi and W 2 are polyno¬ 
mial functions with respect to InlXj^j: for i = 1 or z = 2, Wi{t, X^, X 2 ) = '^1=o't^gi,s{X^jX^)- The 
functions ^ belong to L^{B) and are such that ||(/i,sel+l/^||i 2 (g) < +00 (they are exponentially decaying). 
Therefore, the third point of Proposition 4.2 ensures that i?m, p.bl belongs to SOjj ^ (12+) for /3 < 7 — (Am — p) + |, 
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and consequently for any /? < 2 (Am p\ 


Finally, collecting (D.22)-(D.26) and (D.28), we see that 


AU, 


Tn,p,+ — 2Vx(i? )'V(x 


macro,+^m , macro 


+X-(-^]^)^m,p.BL)+Ax(f?''')(Xi 

+ Rm,p,BLiXi 




n^m. n . macro +X-{X+)U^,p ,bl) 
m,p,macro(^l : X2 )• 


m,p, macro 

R 


which belongs to for for any /? < 2 — (Am — p)- 


D.3 Proof of Proposition 4.5 

The first part of the proof consists in showing that Problem (4.1) has a unique solution u G 03^, ..|,(fl~*’) for any 
/?' < 1 /2. As usual, the proof of existence relies on a variational existence result, which is given by Proposition 4.1 
in the present case. We first verify that the hypothesis on / and g of Proposition (4.5) guarantee that we can apply 
Proposition 4.1. Since A < 7 < 1, > 1, we have 


Vi + {R+¥f 


L2(a+) 


< C\\{l + R+)p^+^ 


L2(n+)- 


Therefore, since /3 > 3,^ < 7 < 1 and /3 — 7 — 1 > 1, + ( 3 ?+)^/ belongs to L^(n+): 




L2(n+) 


< c||(i + < c\\f\u 


(D.30) 


Then, g G iT^^^(rhoie) being compactly supported. Proposition 4.1 ensures that Problem (4.1) has a unique solution 
u G Q3(n'''). Next, we check that the solution u belongs to .y_i(A2+) for any /?' < 1/2: since 76 ( 5 ,1), 

p( 7 -i)-i+p+< 5 p,o < 1 forp = 0 orp = 1 , which means in particular that 


(n+) — ^ 

p—1,7 — 1 ' ' 


{l + R+f 


-7-1, 


L2(a+) 


+ 


{l + R+)^’-'^\/v 


L2(a+) 


For any /3' < 1/2, the right-hand side of the previous equality is controlled by ||'i'||y(Q+p which confirms that 
u G 23^/_i ..y_i(A2+) for any j3' < 1/2. In fact, using a weighted elliptic regularity argument (Proposition 4.1 
in [32]), taking into account that Fhoie is smooth and that g G iT^A^(rhoie) is compactly supported, we deduce that 
u actually belongs to Q3|, .^(12“'') for /3' < 1/2. Moreover the following estimate holds: for any f3' < 1/2, there 
exists a constant C such that 


(a+) 
0 ' 



(D.31) 


It should be noted that the assumption on the smoothness of Fhoie is obviously required for having v G Q3|/ ,^(12+). 
Problem (4.1) being linear, the uniqueness of u S Q3|/ ,^(12+) turns out to be a direct consequence of Proposi¬ 
tion 4.1. 


The second part of the proof consists in proving the asymptotic formula (4.13). It relies on a finite number of 
applications of Lemma 4.6 in order to exhibit the asymptotic behavior of u as i?+ tends to -l-oo. We start with 
applying Lemma 4.6 picking /3i = § < 5 (the choice of 2/5 is not important) and ^2 = /3i + 1 — e, e > 0. Then, 
for e sufficiently small, /3^ is admissible and the sum in (4.15) is empty since there is no integer value of k such 
that 

-?<1-/32<?/c<1-A = | 

As a result, choosing = /3^ — e = | — 2£, we deduce that u is actually in 23^ i_ 2 e(^^) any e > 0 and 
for 1 — 7 sufficiently small. In fact, applying a second time Lemma 4.6 with (3^ = 3 — 2e and an admisible 

/3^ < /3^ -f 2£, we see that u is in 23| (A2+) for 1 — 7 sufficiently small. Of course, the following estimate holds 

5 '7 


Iffl: 


7/5,7 


(n+) 


< C 





ffl/" (Fhoie) 


(D.32) 


Next, let /3^ = I and /3^ = /3^ -f 1 — £ = ^ — e. For e small enough /3^ is admissible, and we have 

<1-/32 <-?< 1-/32. 

Thanks to Lemma 4.6, there exist a constant and a remainder u belonging to 03^^ 7 (^ 2 "*") for any /3° < 0^ 

(providing that 1 — 7 sufficiently small) such that 


u = -I- u with 


|if_l(u)| -f ||u| 




0^ ,■ 


(n+) 


< C 




(n+) 


bll 


(Fhoie) 


49 





















Again, applying a second time Lemma 4.6 with (3^ = ^ — e and = ^ + e for e small enough, we see that u € 
® 12 / 5 . 7 (^’'^) (the sum in (4.15) is empty in this case), with (0+) < ^ {\\f\\m° ^{n+) + ll5llrri/2(fMc))' 

Now, we shall pursue the procedure, working only on the remainder. In fact, we shall consider the slightly different 
remainder 

ui = u — ^-i{u)U-i,p^+, 

where the positive integer p > 3 is chosen sufficiently large in order to ensure that Aui belongs to 
(Lemma 4.3 guarantees the existence of such an integer p (take for instance the smallest integer p such that p > 
P — 2 — |)). Moreover, since 

Ui = ft + Af_i(M)(Z^_1^3_+ — + 

ui belongs to 23^2/5 7 (^^) ~ 7 sufficiently small. We naturally have 

Note also that the restriction gi of 9„tti to Lhoie is compactly supported and belongs to i7^^^(rhoie) since Lhoie is 
smooth. Then, we shall apply again Lemma 4.6 with ^ and P"^ = ^ — e. We remark that, for e sufficiently 
small, /3^ is admissible and, for k = —2, and k = —3, 

1-/3^ <1-/3^. 

Then, Lemma 4.6 ensures that there are two constants Af_ 2 (w) and such that 

Ui = if_2('u)Z^_2,3,+ + if_3(u)W-3,3,+ + “i, Ui € for any /3° < /3^. 

Moreover, 

\-^-2{u)\ + \^-3{u)\ + ||u2||q 32^ (f2+) < ^ (ll/ll®0 ^(^+) H (pMe)) ' 

Applying a second time Lemma 4.6 with P^ = ^ — e and /3^ = ^ + e for e small enough, we see that ui is 
indeed in 93^7/5 and satisfies 

ll^2|la,2^^^^^(fj+) < C (||/|l2jo^^(n+) + ll5llHi/2(rM.)) ■ 

Then, we introduce the new remainder 


U2=Ui- ^-2{u)U-2,p,+ + ^-3{u)U-3^p^+, 

where p > 3 is chosen sufficiently large in order to ensure that Am 2 belongs to 93^ Of course, U 2 belongs 

to 93^7^5 Again, the restriction p 2 of dnU 2 to Lhoie is compactly supported and belongs to iT^^^(rhoie) 

since Lhoie is smooth. 

Repeating the procedure several times (applying Lemma 4.6 at most (2[/3J + 1) times) provides the asymptotic 
formula (4.13) and the associated estimate (4.14). 

E Technical results associated with the matching procedure: proof of 
Lemma 5.1 

We first prove (5.3) by induction. The base case g = 0 is obvious because Wm,o,± is independent of lni?+ and 
Cm, 0,0 = 1 (cf- Remark 5.2). 


Before we prove the inductive step, it is interesting to consider the cases r = 1 and r = 2 separately. For r = 1, 
we remark that Wm,i,+ (6*^, Ini?"*") does not depend on Ini?+. Indeed, Am — 1 = — 1 does not belong to A 

(Am — 1 is not a singular exponent). Then, Lemma C.l ensures that , lni?+) = Wm,i,+ {S'^)- It follows 

that u;m,i,+ (0’^, lni?+ + InJ) = lni?+). By contrast. For r = 2 and m i, Xm — 2. 

belongs to A. More precisely Am — 2 = Am- 3 - Then, we know that 

Wm,2,+ {0"^ MR^) = Wm,2fi,+ {^^) + Wm, 2 . 1 ,+(6'’^) In i?+, with Wm,2,l,+ {6^) = C Wm-3.0.+(^'^)■ 

and C = ^ fg~ Wm, 2 . 1 .+ (('+), Wm- 3 .o,+ (('^)c(('+- Therefore, 

Wm.2,+ (6'’^,lni?^ +ln(5) = Wm.2,+ (6'^, lni?+) + C'm.2.1,+Wm-3,0.+ (6''^) ln5, 
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which proves formula (5.3) for r = 2. The previous argument may then be repeated inductively. The key point to 
figure out is that for any g G N such that g 7 ^ m/3, Am — 2q is always a singular exponent, which corresponds to 
Am—3(?- 


Let us start the inductive step. Let r G N, and assume that formula (5.3) is valid up to order r — 1. Let us denote 
by and the left and right-hand sides of (5.3), i. e., 

+ln(5), 

L/2J / k \ 

Wm-2.k,r-2k,+ {0^ MR^) ( C'm,2fc,t,± (In^)" j . 




ki=0 


wj^HS ''^RHS ^ polynomial dependance with respect to In <5 and In i?+. We shall see that = 

(^+)A„- 

^''JJrhs satisfy the same problem and apply Lemma C.l to conclude. 

First, we can rewrite v^j^rs as 

vlRs{0+,lnR+) = 6^-^- (0+, In (<^^(i?+))) with /(i?+) = R+S. 

Then, in view of the formula {R'^dR+) {g{ip^{R'^))} = ip^{R'^)g'{ip^{R'^)) and reminding that the function 
(p^™'~'^Wm,r,+ i0^ is harmonic in and /C+’^, we have 

Avrhs = Avf f 7 c; = 0 in U 


‘LHS 


Then, we evaluate separately the jump values of and w^Rg. 

r—1 1 — 1 

Kffs] =J2'^r-p9ln,P,r-p,+ i^^R^ + ln<5) + J2'^r-p9l,p,r-p,+ i’^^ >5) 

P—0 P—0 

But, it is verified that 

f)j 

iR+)^--^-^gRj^+(\nR++\n6) = {-ly lni?+ + ln<5))) . 

Then, for p < r — 1, using the induction hypothesis and the fact that Am — P = Am-Sg — (p — 2g), we have 

b/2j q 

9n,p,r—p,+ {^^R T In J) ^ ^ ^ ^ (i^ ffm—3Q.iO—2q,r—p.+ (i^ R )’ 

q=0 i=0 

and, analogously. 


(E.l) 


Lp/2J q 

p",p,,._p,+ (lni?++ln(5) = XI Xl(^^'^)*‘^™. 29 A.+ ffm-39,p-2g.r-p,-r(lni?+) ip<r-l). (E.2) 

q—0 i—0 

Therefore, 

r-l Lp/2J q 

Kr/s]=X J2 (22/_pg^_3,,p_2g,..-p.+ (lni?+) + T>/_pPm-39.p-2«,r-p.+ (lni?+)) • 

p^O g—0 z—0 

(E.3) 

On the other hand, 

Lr/2J q 

= 'y ' y Cm,2q,i,+ [Wm-3qp —2p,-|-] 

q—0 i—0 

[r/2j q r-2q-l 

^ ^ ^ ^ ^ ^ (}^^y^'m,2q,i,-\-(y'^r—2q—p9m—3q,p,r—2q—p,-\-{^^^ ) ^r—2g'—p^m—3g,p,r—2g'—p,+)) 
q—0 i—0 p—0 

Then, using the change of index p p -\- 2q, and interchanging the sums over p and q, we see that [wrhs] = 
[wlhs]- Similar arguments yield [dg+WRRs] = [dg+WLHs]- Einally 

A{vlhs-vrhs) =0, [d 0 +{wLHS - wrhs)] = 0, and [wlhs - wrhs] = 0- 

To conclude we consider separately the case r odd or r = ^ from the case and r even and r 7 ^ ^. If r is odd 
or r = Am — r is not a singular exponent. Then, Lemma C.l ensures that vj^Rg = ''^rrs- if ^ is even and 
r y then Lemma C.l guarantees that v^Rg — Vrrs is proportional to tUm- 3 r/ 2 ,o,+- in other words, writing 


L/2J 


L/2J 


inf?’*’) — X tiffs'.! (i’a and Ini?''") — X ’’^rhsA^^) RT > 


i=0 


i=0 
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we know that for z 7 ^ 0, and WpppgQ = w^p^jpgQ + C{6)wm-r/ 2 ,o,+- We shall see that C((5) 

vanishes. Indeed, 


/ wiHgoie+)w^_3r/2,O,+ i^'^)d0~^ = / Wm,rA^'^)Wm-3r/2,oi^~^)d&~^ 

Jo Jo 

Q_ 

i=0 

On the other hand, using the orthogonality condition (C.15) for fc < r/2, fc 7 ^ to/ 3, the fact that Cm,2k,i,± = 0 
for any integer i such that 0<i<k if k = m/3 (See Remark 5.2), we see that 

/ WgHsi^'^)'^ni-3r/2,0,+ {^^)dd^^'^'^C!m,2k,i,+ {^n6y / Wjn-3k,T-2k,0,+ {S''')Wm-3r/2,0,+ {^'^)dO'^ 

k=0 i=0 "^0 

lr/2\ 3^ „ r/2 

= <^m,r,^M^nSy / \w^-3r/2,oA^A\^de+ = ^ ^ (In<))\ (E.5) 

i=0 do 

As a consequence Wpppg g = w^n^g g and the proof of formula (5.3) is complete. 

Finally, Formula (5.4) then results from (F.1)-(F.2) (which finally hold for any {n,p, r) G N^), and (3.32): 


p„,.,+ (ln(5|X+|),X+) = ^ (gl,,_p,p_+(ln(J|X+|))W^;(X+) +g(;,,_p,p,+ (ln(5|X+|))VF;(X+)) 

p—0 

r [{r-p)/2] q 

= E E (ff™-3,..-p-2,.p.+ %‘(X+) + 5:/-3 ,..-p-2,.p.+ %"(X+)) 

P—0 q—O i—0 

Yr/2\ / q \ /r-2q 

— E! I ^m,2q,i,+ j ( E! ffTn- 3 o.r-p- 2 o.p. + W/, (X^) + gm-3o.r-2o-p.p. +W/, (X^) 

q=0 \i=0 / V p=0 


= E P™-39..-29.+ (ln|X+|,X+) E(ln^)*C'm.29.*,+ . 


F Technical results for the justification of the asymptotic expansion 

El Evaluation of the modeling error: proof of Lemma 6.7 

On the one hand, a direct computation shows that 

E <5*”+^A<,= E <5i"+«-2(Ax^^,, + 20,,9x,^^,, + C^n.,-2)(xl,E+7^'^„,BJx) 


(n,9)er>jvo 


(n,g)eDNo 


where 


L^J a „+2 


T^%0,blA) = E + E E a„ = [No - (F.2) 

n—0 n—0 q—Oiri-A^ 

On the other hand, writing a Taylor expansion of uf^ ^ for ±X 2 > 0, 

°‘^-q f)PyS p pX2 f)a„-q+l.,S ( _ p\a„-q 

.4,fe,x,)= g ^(x.,0*)^ + / dt. (E3) 

we see that 


y: «i"+.|A.x(^)K,,(x)= y; jin+,-»y3^^^(x,,o±)(|Ax,x±W|Si(i) 

(n,9)eDjVo (n,q)&DNg ± P=0 2 \ / 


+ ^^o.taylor(x) (F.4) 
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where 




'Nq ,taylor 


w = E E 




Q^a„-q+l 


-ixi,t)- 


(an - qY- 


dt. 


± {n,q)^DMo 

Subtracting equalities (F.l) and (F.4), taking into account the periodic corrector equations (1.15), we obtain 

fmod,2 = -(1 - xtW - X-(x))ln^_^, (^Aro.taylor(x) + 'R-No,Bl{^)) ■ 


(F.5) 


(F.6) 


It remains to estimate tayior(^) blY^) ^mod- This should be done precisely because the two terms 

are singular close to the two corners. For 7?.^^ tayior(^)’ prove the following estimate; 

Lemma F.l. For any e > 0 small enough, there are two positive constant C > 0 and Jq > 0 such that, for any 

S < So, 

6 


(1 - X+W - X-(x))ln^ T^No,laylori^) 






v{S)) 


(F.7) 


Proof. We first note that, for (5 sufficiently small, the functions Ifjs [A, x±(^)] (^ii 


are supported in the union of the bands By and B 2 where 

By = {-L + p{5)/2, L - 77((5)/2) x (5, 25) B 2 = {-L + p(J)/2, L - p{5)l2) x (-2^, -5). 


{a„-q)\ 


In what follows, we shall obtain an estimate for 
similar analysis can be carried out for i? 2 . 


(1 - Xi(x) - X-(x))ln;5^^7^^„,taylor(x) 


L^Bi 


but a strictly 


In order to separate the singular (close to the corners) and the regular (far from the comer) behaviors of the 
macroscopic terms, we localize the error, namely 

(1 - xtW - X-(x))1^^^^7^^„,,aylo.■(x)||^^^^^^ < C |^|| 7 ^^o,taylor(x)||^,( 50 ) + E 11 .taylor W 11 (b± ) ^ , 


where 


Bt = i-L + p{5)/2, -L/2) X (0, 25), Bf = (L/2, L - p{5)/2) x (0, 25), 73? = (-L/2, L/2) x (0, 25). 


We start with the analysis of 


7^ 


A^o daylor (x) 


L2(B?) 

tation, we see that there exists a constant C > 0 (independent of ) such that 


, where the functions ^ do not blow up. By a direct compu- 


[A,x±(y)] 




Qiar,-q+l 


-{xi,t)- 


(a„ - q)\ 


dxidx2 < C5^°‘" 


(F.8) 

But, a standard elliptic regularity estimates (see e. g. [10]) shows that is smooth in By. Reminding that ^ 
has a polynomial dependence with respect to In 5, it is verified that, for any fc S N, for any e > 0, there exists a 
constant C such that ||u^_q||_H'fc(B 0 ) < C5~‘^. Summing (F.8) over (n, g) € Djvq, using the fact that |n—[|nj > 0, 
we see that 


7^' 


r(x) 


Afo.taylorV^yil 




(F.9) 


We turn to the estimation of 

2n 


'^Aro,taylor(x) 


.We remind that the restriction ut „ to FIt belongs to Vf o (F2t) 

L2(B + ) 

where I3n,q = 1 + ^ + g + e (£ > 0). Using a weighted elliptic regularity argument (see Corollary 6.3.3 in [26]), 
q belongs to ^^^(Ut), for any .( S N, and, as a direct consequence is in ^_ 2 +q(^t)- The 

estimation (F.8) (replacing B^ with 73+) remains valid unless and 

are not uniformly bounded anymore. Nevertheless, since ^ + e + «„ > TVq, 


[ \d)^r‘‘^^^iq\^dx= [ |a“;-«+iw^/(r+)2(¥+^+“A(^+)-2(^+e+a„)rfx< C'??(5)-2(^+^+“-). 
Jb+ ’ JBt 


Introducing (F.IO) into (F.8) and summing over (n, q) & yields 


(F.IO) 


7^' 


Nq ,taylor 




(u,q)(yD^^ 


+E+a„ 


(F.ll) 
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Since 6 < ri{S) and that ^+e + a„=NQ + ^ — > 0, we end up with 




rW 






ri{S) 


No+e 


(F.12) 


□ 


Naturally a similar estimation holds for . Collecting the results of (F.9), (FI2) gives the desired result. 

The analysis of the (x) leads to the following lemma: 

Lemma F.2. For any e > 0 small enough, there are two positive constant C > 0 and 5o > 0 such that, for any 

S < (5o, 

/ § \ ^0+3/2 


(1 - X+(x) - X-(x)) lo^ 




s :-£-2 


(F.13) 


Proof. As in the proof of Lemma F. 1, we shall decompose the error into three terms, which corresponds to evaluate 
the error in three distinct regions of (where 11 ^ ^ are either regular or singular): 

(1 - X+(x) - X-(x)) ^ W||L"(O„od.0) + \\T^n„,bl{^)\\lhqL, ^) 


where 


L 


^mod.O = 1 (a;i,a;2) € “TT < l^^ll < 


^Ld,± = <1 (XI,X2) G <±iL- Xi) < ^ 


We start with the estimation of BL(^)llL 2 (n^ ^ q)- view of the definition (2.30), 


d^,dx,Uiq{xi,^) = ^9S+^(<,g_p(a;i,0))r5xilFp*(^)+^aP^(5^2M^_,_p(a;i,0))r9xiVFp"(^). (F.14) 


p=0 


P=1 


so that 




P—0 


p^l 

where Fmod.o = {(a:i, 3 : 2 ) G F, \xi \ < ^}. In function ^ and its traces over q n F are smooth 

(they have a polynomial dependance with respect to In 6). Consequently, for any A: G N, for any e > 0, there exists 
a constant C > 0 such that 

l|5S;^^(Mi,5-p(a^G0))r||L~(r„„,.o) < and 0))r||L=o(r„„,„) < (7(5-®. (F.16) 

Then, using the periodicity of the boundary layer terms, we obtain 


|- 2L-2 ^ 


/ |5x,VF;(-)pdx< ^ 6^dx,WXLHB)<C6, 

'^mod,0 £—0 


Similarly, 


|lL;(J)| 2 dx <52 ^ \\w;\\Ihb) < cs. 


(F.17) 


(F.18) 


^^mod,0 "" £—0 

Collecting the results of (F. 14)-(F.15)-(F.16)-(F.17)-(F. 18), we obtain 

A similar analysis leads to 

Then, summing over n and q in (F.2) gives 


(F.19) 


We can now evaluate bl(x)||l 2 (q 5 (the evaluation of BL(x)||B 2 (f 26 ^ ) being similar). Here 

again, the inequality (F.14) (replacing Fmod.o with rniod,+ = {(a:,i ,^2) G r,p(5) < xi < f } and H^^d.o with 


54 



^mod +) valid, but the norm of the tangential traces of the macroscopic helds are not uniformly bounded any- 


^Q+ly5 


'n,q belonging to both and 

+„+i(^b) (/3 n,9 = 1+^+g + e). Consequently 9“^ and belong to 


more. However, for any positive integer a > 0, d'^^dx^u^ q and 9““ 

T/2 

which means in particular that (see Lemma 6.1.2 in [26]) that || |a;i — L|^".9+“+i/2£)“^ {dx^^n g)r||L°°(r„„d +) and 
II |a:i — L|/5 ".‘!+“+i/2(^“^ (5 a:2'“n.q)r||L“(r„od +) uniformly bounded. As a consequence. 


Finally, 

Similarly, 




< C5i/2r7(<5)-(5+^+T+«). 

Then, summing over n and q in (F.2) gives 


/ S \ -^0+3/2 

ll^k.sLWIlL2(ni,,+) < C* 


(F.20) 

Collecting (F.19) and (F.20) finishes the proof. □ 

Finally, the estimate of the modeling etTor of Lemma 6.7 results from (F.6), Lemma F. 1 and Lemma F.2. 


F.2 Evaluation of the matching error 

F.2.1 Proof of Lemma 6.9 

Since TZmRcm,n,q,s G V" 2 ^^(Ht) n C2^^(f2B) for any /? > 1 — ^(fe("’9^^o)+i) ^ using the fact that k{n,q,No) > 
|(A'o — q) — n — ^ and the fact that r+ < 2q{6) on one can verify that, for any £ > 0, there is a positive 

constant C > 0 such that 


|7^ 


macro,n,(3,(5|li2(Q+ ) < V{^) 

' match 


No-in-q-e+i 


|V7^ 


•macro,n,g,(5 | 


< v{S) 


No-in-q-e+l 


(F.21) 


As a consequence, summing over (n, q) G D^o remarking that 5 < rj{S) and Vxmacro.+ (x'''/^) < 5 ^ completes 
the proof. 


F.2.2 Proof of Lemma 6.10 

The function TZgj^ is supported in the domain —L + 5 < Xi < L — S. We shall evaluate it in the domain 

^match.BL = ^match ^ 2 : 2 ) G K ,Xi — L < —(5}. 

that we separate into two parts ^ and 2 defined as follows; 

^match.BL.i = ^match ^ {(a^i, a;2) G ^ < a^i - L < -<5} 

^match,BL.2 = ^inatch ^ {[Xl,X2) G —2p((5) < Xi — L < —} 


We shall now estimate TlBL,n,q,s (or, more precisely we study each term of the form {wn,q,j{xi,0))Wn,q,j) in 
^inatch BL 1 ^match BL 2 ' ’^^e domain the profile functions Wn,q,j are exponentially decaying 

(since |a: 2 | > ^q{5) on this domain). As a consequence, for any N G'H there exists a constant C > 0 such that 


|7^' 


BL.tVollL2(n+,^ 


|V7^' 


BL.Afol 




(F.22) 


Next, we remark that for any/3 > 1 —^ is bounded ((r(;„_qj(a;i, 0)) G 
^ 2 ^|s+i(^))- ^ consequence,for any £ > 0, 

' match •' 

Then, using the previous estimation and the periodicity of the profile function Wn,q,j, we have 


Wn,q,j{Xi,Q))Wn,q,j\\B'^(Q,+ ) ^ C’|| (tUn.ij.j , 0)) || /p+ j||bF, 

^ match, is L. 2'• match-' 


'".9jilB2(n+ 

< C'(55p(<5) 


match, BL .2 ■' 

No+h-q-i 


(F.23) 
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Analogously, for any/3 > 2 fc(ri,g Afo)+i ^ ||| 3 .^_ 7^|^+2 (ifn,g,j(a;i, 0))||^oo(p+ 1 is bounded. As a consequence, 

' match/ 

for any £ > 0 , 

IIV ((u;„,,,,(xi,0))W^„.,.,) 11^2(0+,(F.24) 
Summing up over j G (0, K) and (n, q) G Dn^, noting that x- C^V^ ) = 1 on 2 ’ obtain. 


|7^ 




Bi.WollL^(a+ 2 ) -^ ii'"''SAiVollL2(n_, ,, 


|V7^f 


< (F.25) 


Collecting (F22) and (F.25), and carrying out an entirely similar analysis for \\T^%l 1 ^ 2 (^qq+ pp'*) l^^ds to 
(6.36). 


F.2.3 Proof of Lemma 6.11 

As previously, we first investigate separately 'R-MF,n,q- We consider the scaled matching domains 




t){5) 


MS) 


match 


= {X+ e ^ < i?+ < 2 ^}, rLeh = {X+ e ^ < i?+ < 2 ^}. 


viS) 


MS)- 


Then, making the change of scale X+ = x+/(5, we have 

l!'7^tVF,n,(j||p2fo+ ) ^ S^ f iT^ATFn 5pC?X+ 




< (52||(l + i?+)-('^-'t-l)p-(T'-l)||2 


If 7 — 1 < 0, p < (1 + R^) t). As a result, for any /3 < 1 + 2L“»»^gJ+i ^ ^ g i ^ ^ 1 sufficiently 

small, 

/ (5 ^ ^ " 

\\RNF,n,q\\L2(^^^j < C6 

Taking /3 = | + [Nq — g) — |n — 1 — e, e > 0 (| [ctn^gj > {Nq — g) — |n — 1), we have 


\MNF,n,q\\L2^^+^ < CS 


No—q—4n—2—e 


Summing over (n, q) € PjVo^ we obtain 


|7^ 


AfF,AfollL2(o4 


<C 5i"+«<5 

(n,(})GFjvo 


f — 


-^0 —9— —2 —£ 


< C(5 


77(5) 


Afo-2-e 


The gradient of TZ%p can be estimated in the same way noticing that 

\\'^'R-NF,n,q\\‘l2(n+ ) ^ [ IVx+'7^]VF.n,gPdX+ 

V^^match-^ /0'5 

^ ^match 

< ||(l + i?+)-(''-^V-^'^-'^lli=o(A^^_Jll(l + ^+)''-V-'Vx+7^^F.„.,|li2(a+) 

As for the trace estimate, we note that \\{R^)'^TZNF,n,q\\fq 2 (QS \ is uniformly bounded for any 

''^^match/ ^ 

Consequently, since F^^tch smooth, under the same condition for 7 , ||7^ArFn,q(.R~'')^|| rcxj^f^ i is uniformly 
bounded. Then 

1 1 / (5 \ ^ 1 / (5 

ll'7^iVF.n,q|lL2(af2+^^__nr^) < < C6^ j - 


- 1 . 
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